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PREFACE 


This  report  is  written  in  fulfillment  of  the  requirements  of  U.  S.  Navy 
Contract  NOw  66-032*+c.  The  contract  is  a  "level  of  effort"  type,  and  the 
objectives  are  covered  by  the  following  work  statement,  taken  from  the  contract. 

"Develop  and  verify  an  accurate,  general,  and  rapid  method  of  calculating 
axially  symmetric  and  two-dimensional  turbulent  boundary -layer  flows.  The 
specific  phases  of  the  work  are  as  follows: 

(a)  Semi -empirical  expressions  will  be  developed  for  the  turbulent  trans¬ 
port  properties,  such  as  eddy  viscosity,  for  use  in  solution  of  the  flow 
equations. 

(b)  Solutions  of  the  complete  partial  differential  equations  will  be  ob¬ 
tained  for  incompressible  flow,  covering  problems  of  flow  in  water  and  in  air. 

(c)  Solutions  will  also  be  obtained  for  compressible  flow,  applicable  to 

air. 


(d)  Illustrative  applications  of  the  method  will  be  performed  to  provide 
information  on  velocity  profiles,  boundary-layer  thickness,  skin  friction  and 
heat  transfer. 

(e)  Accuracy  will  be  checked  by  solution  of  a  variety  of  flow  problems 
and  comparison  with  experimental  data." 

Under  this  type  of  contract,  when  the  due  date  arrives,  accomplishments 
are  reported  whether  the  ultimate  goals  have  been  reached  or  not.  In  the 
present  case,  as  might  be  expected,  the  studies  are  not  complete.  The  problem 
of  incompressible  flow  has  been  rather  well  explored,  although  more  remains  to 
be  done.  The  results  have  been  surprisingly  good.  The  equations  governing 
compressible  flow  have  all  been  programmed  and  the  method  is  working,  but  time 
was  available  to  calculate  only  a  few  cases  of  flat -plate  flow,  with  and  without 
heat  transfer.  Much  more  work  remains  to  be  done,  and  the  gratifying  results 
for  incompressible  flow  supply  a  firm  foundation  for  continuation,  which  should 
follow.  The  present  method  has  been  programmed  on  the  IBM  709**  unler  the  num¬ 
ber  57EB.  The  program  can  be  obtained  by  qualified  requesters  from  "Commander, 
Naval  Ordnance  Laboratory,  White  Oak,  (Code  330 ),  via  Commander,  Naval  Ordnance 
Systems  Command  (Code  0RD-035)"* 

The  authors  .nd  their  company  wish  to  express  gratitude  for  the  support 
supplied  by  the  U.  S.  Navy.  Without  it,  it  is  unlikely  the  work  would  ever  have 
been  accomplished. 
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1.0  SUMMARY 

This  report  presents  a  numerical  solution  of  turbulent  boundary-layer 
equations  for  both  compressible  and  incompressible  flows.  An  eddy  viscosity 
concept  is  used  to  eliminate  the  Reynolds  shear-stress  term,  and  an  eddy- 
conductivity  concept  is  used  to  eliminate  the  time  mean  of  the  product  of 
fluctuating  velocity  and  temperature.  The  turbulent  boundary  layer  is  regarded 
as  a  composite  layer  consisting  of  inner  and  outer  regions,  and  a  separate 
expression  for  eddy  viscosity  is  used  in  each  region.  The  ratio  of  eddy- 
viscosity  to  eddy  conductivity  is  assumed  to  be  constant.  An  implicit  finite - 
difference  method  is  used  in  the  solution  of  both  momentum  and  energy  equations 
after  they  are  linearized. 

A  variety  of  flows  have  been  computed  by  this  method,  and  comparisons 
with  experimental  data  and  various  correlations  have  been  very  encouraging. 

The  results  described  in  this  report  do  not  represent  a  finished  development 
but  onxy  what  has  already  been  accomplished  by  using  one  particular  formula¬ 
tion  of  eddy  viscosity  and  constant  turbulent  Prandtl  number. 
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4.0  PRINCIBAL  NOTATION 

local  akin-friction  coefficient,  eq. (6.113) 

average  skin-friction  coefficient 
specific  heat  at  constant  pressure 

local  shear-stress  coefficient  for  laminar  flow,  eq. (6.114) 

local  shear-stress  coefficient  for  turbulent  flow,  ^q. (6.115) 

P  P  ,  viscosity-density  parameter 
pe*e  * 

equilibrium  boundary-layer  parameter,  — 

Tw 

dimensionless  stream  function,  eq. (6.46) 

dimensionless  total -enthalpy  ratio,  eq.  (6.54),  where  applicable 
defect-shape  factor,  eq. (6.122) 
specific  enthalpy 

total  enthalpy,  h+  -  u  or  shape  factor,  eq.  (6.121),  where  applicable 

variable -grid -system  parameter 

mixing  length 

reference  body  length 

Mach  number 

pressure 

Prandtl  number 

local  heat-transfer  rate  per  unit  area 
radial  distance  from  axis  of  revolution 
radius  of  body  of  revolution 
Reynolds  number,  ugx/ve 
Reynolds  numbtf.  u  9/v 

c  C 

Stanton  num’ber,  eq.(6.1l8) 
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transverse  -curvature  terra,  eq.  (6.42),  where  applicable 

absolute  temperature 

x-component  of  velocity 

friction  velocity,  VvT 

y-component  of  velocity 

distance  along  surface  measured  from  leading  edge  or  from 
stagnation  point 
distance  normal  to  x 

angle  between  normal  to  the  surface  y  and  the  radius  r  ,  fig.  4 

dimensionless  velocity-gradient  term,  eq. (6.53). 

intermittency  factor,  eq. (6.26),  or  convergence  criterion, 

where  applicable 

boundary-layer  thickness 

defect-displacement  thickness,  eq.  (6.123) 

eddy  viscosity 

transformed  y-coordinate 

momentum  thickness,  eq.  (6.111) 

thermal  conductivity 

dynamic  viscosity 

kinematic  viscosity 

transformed  x-coordinate 

density 

shear  stress 

perturbation  quantity,  f  —  f 
stream  function 
vorticity 


I 
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SUBSCRIPTS 

c  evaluated  at  the  switching  point  of  the  boundary  layer 

e  evaluated  at  outer  edge  of  boundary  layer 

w  evaluated  at  wall 

oo  evaluated  at  free -stream  or  reference  conditions 
Primes  on  f  denote  differentiation  with  respect  to  t) 
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5.0  INTRODUCTION 

The  boundary-layer  concept,  first  introduced  in  190^  by  Prandtl,  divides 
the  flow  past  a  body  into  two  regions:  an  inviscid  region,  governed  by  the 
Euler  equations  of  motion,  and  a  thin  viscous  region  in  the  neighborhood  of  the 
body,  governed  by  the  boundary-layer  equations.  For  laminar  flow,  the  existence 
of  a  known  relationship  between  the  shear  stress  and  the  velocity  gradient 
completes  a  set  of  partial  differential  equations,  and  exact  solution  of  the 
boundary-layer  equations  is  mathematically  possible.  Highly  accurate  solutions 
exist  for  some  simple  flows,  such  as  similar  flows,  which  are  especially  impor¬ 
tant.  With  the  advent  of  high-speed  computers,  quite  satisfactory  results  for 
a  variety  of  general  flows  have  been  obtained. 

For  turbulent  flows,  on  the  other  hand,  because  of  the  limited  understand¬ 
ing  of  the  turbulent  process,  the  exact  solutions  of  the  boundary-layer  equa¬ 
tions  are  not  possible.  The  usual  boundary-layer  equations  for  such  flows  con¬ 
tain  a  term  involving  the  time  mean  of  the  product  of  two  fluctuating  velocities, 
which  is  known  as  the  turbulent  shear  stress,  and  a  term  involving  the  time 
mean  of  the  product  of  a  fluctuating  velocity  and  a  fluctuating  temperature. 

At  the  present,  these  terms  have  not  been  rigorously  related  to  the  mean  velocity 
and  mean  temperature  distributions.  Thus,  exact  solutions  of  the  boundary-layer 
equations  for  turbulent  flows  are  not  possible.  In  order  to  proceed  at  all, 
the  solutions  must  depend  on  some  empirical  information.  Even  then  the  solution 
of  boundary-layer  equations  is  not  easy.  For  this  reason,  most  of  the  work  on 
turbulent  flows  has  been  centered  on  empirical  correlations  together  with  inte¬ 
gral  methods.  In  general,  the  approaches  followed  in  these  methods  vary  widely. 
In  one  approach,  for  example,  Head*s  method  [1],  the  boundary-layer  parameters 
are  obtained  by  solving  the  momentum  integral  equation  with  two  empirical  ex¬ 
pressions  called  auxiliary  equations.  These  equations  consist  of  an  expression 
for  local  skin-friction  coefficient  (c^.)  and  an  expression  for  shape  factor  (H). 
In  another  approach,  for  example,  Truckenbrodt's  method  [2],  the  boundary-layer 
parameters  are  obtained  by  solving  both  the  momentum  and  energy  integral  equa¬ 
tions  by  using  an  empirical  expression  for  the  dissipation  integral  and  by 
introducing  further  approximations  in  the  solution.  These  methods  were  recently 
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re  viewed  by  Hiompson  [3]  for  two-dimensional  incompressible  turbulent  flows 
and  were  found  to  give  widely  differing  and  often  inaccurate  results. 

A  more  fundamental  approach  to  the  solution  of  turbulent  boundary  layers 
is  to  regard  the  turbulent  boundary  layer  as  a  composite  layer  and  to  charac¬ 
terize  it  by  inner  and  outer  regions  (see  figure  l).  The  existence  of  two 


Figure  1.-A  turbulent-boundary-layer  velocity  profile. 

regions  is  due  to  the  different  response  to  shear  and  pressure  gradient  by  the 
fluid  near  the  wall.  The  inner  region,  whose  thickness  is  approximately  0.1  to 
0.2  6,  depends  primarily  on  the  wall  shear  stress  and  fluid  viscosity.  The 
mean  velocity  distribution  in  this  region  responds  rapidly  to  "hanges  in  these 
wall  conditions  because  the  eddies  in  this  region  are  very  small.  The  mean 
velocity  distribution  may  be  described  by  the  so-called  "law  of  the  wall": 

^  (yu*A0 

This  relation  was  originally  obtained  by  Prandtl  from  a  mixing-length 
concept  [4],  In  addition,  if  an  expression  for  eddy  viscosity  is  introduced 
in  the  inner  region,  it  can  be  shown  that  eddy  viscosity  in  this  region  varies 
almost  linearly  with  distance. 
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In  the  case  of  a  smooth  wall,  the  inner  region  contains  a  layer,  commonly 
called  the  laminar  sublayer,  adjacent  to  the  wall,  where  the  flow  is  primarily 
viscous  and  the  mean  velocity  increases  linearly  with  distance  from  the  wall. 
The  thickness  of  this  layer  is  of  the  order  of  0.001  to  0.01  6. 


The  outer  region,  on  the  other  hand,  contains  80  to  90  percent  of  the 
boundary  layer  thickness.  The  flow  in  this  region  is  independent  of  the  fluid 
viscosity,  but  is  dependent  on  the  wall  shear  stress,  and  it  is  highly  affected 
by  conditions  in  the  free  stream  such  as  streamwise  pressure  gradient.  Kie 
mean  velocity  distribution  is  conveniently  described  by  the  so-called  "velocity- 
defect  law": 


u  —  u 
e 


=  <p5(y/6) 


The  flow  in  the  outer  region  shows  some  similarity  to  wake  flow.  Near  the 
outer  edge,  it  has  an  intermittent  character.  The  turbulence  is  characterized 
ty  large  eddies.  The  response  of  the  mean  velocity  distribution  to  changes  in 
its  determining  conditions  is  much  slower  than  that  of  the  inner  region.  In 
addition,  an  eddy  viscosity,  if  introduced,  shows  a  nearly  constant  value  across 
the  region.  For  example,  as  suggested  by  Clauser  [5], the  eddy  viscosity  for 
the  so-called  "equilibrium"  boundary  layers  is 

* 

c  =  k  p  u  6 

d  e 

where  k^  was  empirically  determined  to  be  0.018. 

The  approach  in  whi^h  the  turbulent  boundary  layer  is  regarded  as  a  com¬ 
posite  layer  consisting  of  inner  and  outer  regions  was  followed  in[6]  and  [7] 
for  incompressible  flows.  In  both  references,  the  Reynolds  shear-stress  term 
was  eliminated  through  the  use  of  an  eddy-viscosity  concept.  The  main  differ¬ 
ence  between  the  two  approaches  is  the  expression  used  for  eddy  viscosity  in 
each  region.  Another  difference  is  the  transformation  used  to  stretch  the  co¬ 
ordinate  normal  to  the  flow  direction  to  reduce  the  variation  of  the  boundary- 
layer  thickness  and  to  remove  the  singularity  at  the  leading  edge  or  at  the 
stagnation  point.  A  third  difference  is  the  method  used  to  solve  the  boundary- 


l 
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layer  equations.  In  reference  [7]  the  momentum  equation  was  solved  in  its  non¬ 
linear  form  by  an  integration  technique;  in  reference  [6]  the  momentum  equation 
was  solved  in  a  linearized  form. 

The  approach  used  in  this  report  is  the  one  used  in  reference  [7].  Again 
the  Reynolds  shear-stress  term  in  the  momentum  equation  is  eliminated  through  the 
use  of  an  eddy-viscosity  concept,  and  a  separate  expression  "or  eddy  viscosity 
is  used  in  each  region  (see  figure  2).  However,  this  cime  the  expressions  are 


Figurt  2.- Eddy-viscosity  distribution  across  a  boundary  layer. 

slightly  modified,  tc  account  for  the  compressibility  effect.  In  addition, 
the  time  mean  of  the  product  of  a  fluctuating  velocity  and  a  fluctuating 
temperature  in  the  energy  equation  is  eliminated  through  the  use  of  an  eddy- 
conductivity  concept  and  is  introduced  into  the  energy  equation  through  the 
definition  of  turbulent  Prandtl  number: 

c  e 

As  an  initial  step,  the  turbulent  Prandtl  '.umber  is  assumed  to  be  constant. 
Note  that  the  present  framework  is  general  and  tliat  it  can  handle  widely  vary¬ 
ing  eddy-viscosity  and  turbulent  Prandtl  number  formulations. 


The  method  used  to  solve  the  boundary-layer  equations  here  is  different 
from  the  one  used  in  previous  studies  [7  through  13].  Again  the  streamvl  «?e 
derivatives  in  both  momentum  and  energy  equations  are  replaced  by  finite 
differences.  However,  unlike  the  previous  studies,  in  this  study  the  momentum 
equation  is  solved  in  a  linearized  form,  and  the  previous  integration  tech¬ 
nique  is  replaced  by  am  implicit  finite-difference  technique.  A  variety  of 
flows  has  been  computed  by  this  method,  amd  comparisons  with  experimental  data 
and  various  established  correlations  are  made.  The  results  presented  in  this 
report  do  not  represent  a  finished  development,  but  are  only  the  results  ob¬ 
tained  so  far  by  one  particular  eddy-viscosity  and  I'onstant  turbulent  Prandtl 
number  formulation,  in  fact  the  first  one  tried. 
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6.0  DESCRIPTION  OF  METHOD  OF  SOLUTION 


6.1  Equations  of  the  Compressible  Turbulent  Boundary  Layer 

The  governing  equations  describing  the  flow  about  two-dimensional  and 
axisymmetric  bodies  at  high  Reynolds  numbers  and  constant  pressure  within  the 
boundary  layer  are  [11]: 


CONTINUITY 

St +  •\{sr (rkp  u> +  57 <rke  v)}  ‘ 0 

MOMENTUM 


cti  .  ou  ,  „  du 

pst  +  pusr+  p 


=  -  in  + 
dx 


q.  ^  / 


5y 


(6.1) 


(6.2) 


ENERGY 

dH  . 

dH  . 

dH  1  3  k 

1  u  dH  . 

P  5t 

p  3x  p 

(4 

Vs 

r- 

ii 

vs 

> 

\5T3y  + 

where 

k  =  0  for 

two-dimensional 

flow  and  k 

p(1-'fe)u^}  (6-5) 

=  1  for  axisymmetric  flow. 


The  basic  notation  and  scheme  of  coordinates  are  shown  in  figure  3, 
where  u^  is  a  reference  velocity  and  ug(x)  is  the  velocity  just  outside 
the  boundary  layer.  The  term  He,  which  is  a  constant,  is  the  total  enthalpy 
outside  the  boundary  layer.  Local  enthalpy  outside  the  boundary  layer,  namely, 
h  ,  is  given  by 

'  H  .h  4.2 

e  e  2  e 

The  coordinates  are  a  curvilinear  system  in  which  x  is  distance  along  the 
surface  measured  from  the  stagnation  point  or  leading  edge.  The  dimension  y 
is  measured  normal  to  the  surface.  Within  the  boundary  layer,  the 
velocity  components  in  the  x-  and  y-dire°tions  are  u  and  v,  respectively. 

The  body  radius  is  r  . 

o 


Figure  3.- Boundary  layer  on  a  body  of  revolution.  Coordinate  system. 


In  these  equations,  the  t  ransve rse -curvature  terms,  which  are  of  second 
order,  are  retained  because  of  their  importance  in  preoicting  boundary-layer 
growth  on  long  slender  bodies  such  as  certain  missiles  or  at  the  tail  of  a 
streamlined  body  of  revolution. 


The  equations  (6.1),  (6.2),  and  (6.3)  apply  to  turbulent  as  well  as  to 
laminar  flows,  providing  the  dependent  variables  —  velocity,  density,  and 
enthalpy  —  are  replaced  by  their  inrtantaneous  values  for  turbulent  flow. 

The  procedure  is  due  to  O.Reynolds.  With  instantaneous  values,  (6.1),  (6.2), 
and  (6.3)  become 


*pi  + 

sr  + 


.  i  3  /  k  \  3  /  K 

T  (r  pJui)  *  (r  1 

du 

’ivsr  * 

■i „  aHi ,  a,Ii  1  a 

ci"St  *  +  ViTy  ’  T  ^ 


du. 


pi  ST"  +  p<u< 


du. 


i  dpi 

Vi  W  ’  - 


■pivi  >}  - 0 

i  a 


dx 


i  a  ,  •  k  \  \ 


dH .  .  du .  I 

w +  p(1  _  &  vs^j 


(6.4) 

(6.5) 

(6.6) 
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Let  the  instantaneous  values  be  denoted  by  their  average  and  fluctuating 
values,  as  follows: 


u^  =  u  +  u' 


Pi  =  P  +  P' 


V.  =  V  +  V* 


H.  =  H  +  H* 
1 


=  P  +  P' 


(6.7) 


Introducing  the  first  three  expressions  defined  by  (6.7)  into  (6.4)  and 
averaging  with  respect  to  time  gives 


^  rk(p  u  +  p'u')  +  ^  rk(p  v  +  p'v')  =  0 


(6.8) 


By  time  average  is  meant,  for  example, 


U  =  t  J  Ujfr  +  T )dx 


with  t  large  compared  with  the  time  scale  of  the  turbulent  motions.  For 
the  fluctuating  values,  say  u',  the  average  value,  u*,  is 


^  \  /»’ 


dx  =  0 


The  momentum  equation  for  steady  compressible  turbulent  flow  can  be  obtained  as 
follows.  Multiplying  (6.4)  by  u^  (6.5)  by  r,  and  adding  the  resulting 
equations  gives 


rVi  ST  *  rkpiV57  *  ui  ST  (rk°iui)  +  U1  5?  (rVi) 


b  ,  k 


k  dPi  +  d  ,  k  ^i  s 

=  • r  —  *  57 (r  “ 


(6.9a) 


which  can  be  written  as 


/  k  \  ^  /  k  \ 

r—  (r  p.u.u.  )  +  (r  p.v.u.) 
bx  '  i  1  1'  by  1  i  i' 


k  dpi  a  ,  k  \  , 

■ r  —  +  sr <r  11  w } 


(6.9b) 


Introducing  the  first  four  expressions  defined  by  (6.7)  into  (6.9b),  averaging 
with  respect  to  time,  and  using  (6.8)  gives 
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(p  u  +  p'u*)  ^  +  (p  v  +  p'v’)  ^  ^  ^  ^  rk(p  u'2+  u  p'v'  +  p'u'2  ) 

r 


+  rK{  p  ^  -  “  v'u*  -T  p’u’  "  oW  } 


(6.10) 


The  energy  equation  for  compressible  turbulent  flow  can  be  obtained  in 
a  similar  way.  Multiplying  (6.4)  by  H^,  (6.6)  by  r,  and  adding  the  resulting 
equations  gives 


„  9  ,  k  \  3  /  k  x  ^  k  wi  ^  k  ~‘i 

H1  Si (r  piui>  *  Hi  Si (r  Vi5  +  r  Vi'Si +  r  °iyi  ST 


d  /  k 


d  k 

57r 


dH.  .  du. 

ft  w  *  11(1  ■  *  }  V3F  J  (6-:ua) 


which  can  be  written  as 


Si  (*lc<y>1)(V  +  S  (rVi>(Hi> 


*  t-  f  oH, 

|ir{ftS"+  11(1 


Nsr} 


(6.11b) 


Introducing  the  expressions  defined  by  (6.7)  into  (6.11b),  averaging  with 
respect  to  time,  and  using  (6.8)  gives 

(p  u  +  p'u»)  +  (p  v  +  p'v*)  ^  - - ^  rk(p  u*H'  +  u  p'H*  +  p*u'  H) 

+  "7  ¥  rk{l£  -  P  +  nk  -  ^  -  v  P7^-  P’v’H*} 

(6.12) 

Since  for  flows  at  high  Reynolds  number  the  boundary  layer  is  assumed  to 

be  thin  and  the  terms  such  as  p,  u,  H,  and  x  are  assumed  to  be  of  the  order 

of  1  and  v,  y  of  the  order  of  5,  some  of  the  correlation  terms  involving 

u*,  v',  p*,  and  H*  in  (5.8),  (6.10),  and  (6.12)  can  be  neglected.  The 

double  correlation  terms  such  as  p*u',  u'H,,p,H*,  and  p*v'  are  of  the  order 

of  6  at  most,  and  the  triple  correlation  terms  such  as  p*v,u*  and  p*v,H* 

2 

are  of  the  order  of  6  at  most.  When  these  simplifications  are  introduced 
into  (6.8),  (6.10),  and  (6.12)  and  the  predominant  terms  of  the  same  order  of 


I  1 


- 
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mgnitude  are  retained,  the  resulting  equations  are  of  the  same  form  as  those 
of  laminar  flow  except  for  the  terms  p'v*,  —  pu'^,  —  p  v'u',  and  —  pv'H'.  In 
addition,  if  the  pressure  fluctuations  within  the  boundary  layer  are  considered, 
then  another  term,  namely,  —  ^  p  v*^,  should  be  included  on  the  left-hand 
side  of  (6.10).  For  flow  conditions  away  from  separation,  the  terms  —  p  u'2 
and  —  p  v*2,  known  as  Reynolds  normal  stresses,  are  small  and  will  be  neglected. 
Hence,  with  all  these  simplifications,  the  governing  equations  for  the  compres¬ 
sible  turbulent  boundary  layer  become: 


CONTINUITY 


(rkp  u)  +  ^  rk(p  v  +  p*v*)j  =  0 


(6.13) 


MOMENTUM 


pur  +  (p  v  +  p*v*, 


ENERGY 


dH  /  ■  ■, 

0UST  +  (pv  +  P'y 


7\  dH  1  d  [  k  Ju./.  .  £  Pr  \  dH  .  lx  du  V 1 

=  \*{1  +  11(1 " p?)u 57 ;J 


(6.15) 


where  the  term  —  p  u*v*,  known  as  the  Reynolds  shear  stress,  is  eliminated 
through  the  use  of  Boussinesq's  eddy -viscosity  (e)  concept  and  the  term 
—  p  v*H'  is  eliminated  through  the  use  of  an  eddy-conductivity  (\p)  concept. 


“  P  u»v»  =  e  ^ 
-tuT  ^ 

-  P  V'H'  "  —  57 

P 

C  € 

=  "£r 

T 


(6.16) 


(6.17) 


(6.18) 


If  (  )  denotes  wall,  the  boundary  conditions  to  be  considered  are : 


MOMENTUM 


u(x,  0)  =  0 

v(x,  0)  =  0  or  v(x,  0) 

lim  u(x,  y)  =  u  (x) 
y  ->  00 


vw  (mass  transfer) 


(6.19a) 

(6.19b) 

(6.1Sc) 
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ENERGY 

H(x,  0)  =  Hw  or  (x,  0)  =  (  ^  )w  f<5.20a) 

lim  H(x,  y)  =  H  (x)  (6.20b) 

y  -»  oo 

6.2  Formulation  of  Eddy  Viscosity  and  Turbulent  Prandtl  Number 

In  order  to  solve  the  compre  ssible  turbulent  boundary-layer  equations 
given  in  the  last  section,  it  ?  s  necessary  to  use  expressions  for  eddy  vis¬ 
cosity  and  turbulent  Prandtl  number.  The  eldy -viscosity  formulation  that  will 
be  used  in  this  study  is  the  same  as  the  one  used  for  incompressible  flow  in 
a  previous  study  [7]»  This  formulation  has  worked  well  for  incompressible  flow 
and  hence  it  was  decided  to  extend  it,  with  small  modifications,  to  compressible 
flow.  In  this  formulation,  the  boundary  layer  is  regarded  as  a  composite  layer 
characterized  by  inner  and  outer  regions.  In  the  inner  region,  an  eddy  vis¬ 
cosity  based  on  Prandtl*s  mixing-length  theory  is  used;  in  the  outer  region,  a 
nearly  constant  eddy  viscosity  is  used.  It  is  <  xactly  constant  when  the  flow 
is  incompressible  and  without  heat  transfer.  An  intermittency  factor  is  applied 
to  this  basic  "outer  viscosity". 

6.2,:.  Viscosity  in  the  inner  region. 

In  the  inner  region,  the  eddy  viscosity  is  represented  by  Prandtl's  formula 
based  on  the  mixing-length  theory;  that  is, 

ci =  p  ‘2 1^1  (6-21) 

where  l,  the  mixing  length,  is  given  by  l  =  k^y.  However,  to  account  for 
the  viscous  sublayer  close  to  the  wall,  a  modified  expression  for  l  is  used 
in  (6.21).  This  modification,  suggested  by  Van  Driest  [l^]  and  developed  by 
consideration  of  a  Stokes-type  flow,  is 

/  =  k-jjrfl  -  exp(—  y/A)]  (6.22) 

Substituting  this  expression  for  l  into  (6.21)  gives 
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€i  =  p  kl  y2*'1  ”  exP(“  y/A)]2  |^|  (6.23) 

This  expression,  as  it  stands,  applies  to  incompressible  flows.  The 
quantity  p  is  a  constant,  k^  =  0.4,  and  A  is  a  constant  for  a  given  stream- 
wise  location  in  the  boundary  layer,  defined  as  26v(p/tw)^2.  Equation  (6.23) 
shows  that  as  y  increases,  the  exponential  term  disappears,  leaving  Prandtl's 

form,  namely,  equation  (6.21).  It  also  shows  that,  for  y  -*  0,  e.  should 

U  1 

vary  as  y  ,  The  latter  conclusion  is  in  contrast  to  the  behavior  of  eddy- 

viscosity  expressions  proposed  by  Townsend  [15]  and  Reichardt  [16],  which  show 

that  should  vary  as  y^.  On  the  other  hand,  (6.23)  has  the  same  behavior 

close  to  the  wall  as  Deissler*s  eddy-viscosity  expression  [17];  that  is,  e. 

1+  1 
varies  as  y  as  y  -*  0.  An  analysis  given  in  Appendix.  A  indicates  that  e. 

4  1 

should  vary  as  y  as  y  -*  0. 


Equation  (6.23)  can  also  be  applied  to  compressible  flows  if  p  is  taken 
to  be  a  variable  and  if  the  exponential  term  is  modified  to  account  for  the 
heat  transfer  in  the  sublayer.  A  logical  generalization  is  to  consider  a 
Stokes-type  flow  in  which  the  fluid  has  a  variable  viscosity.  An  analysis 
given  in  Appendix  B  indicates  that  the  eddy-viscosity  formula  for  the  inner 
region  should  now  be 


e 


i 


(6.24) 


l/P  — 

where  A  =  26v  (p^i  )  '  and  v  is  the  moan  value  of  v  obtained  by  averaging 
(6.24)  over  some  arbitrary  distance,  perhaps  the  sublayer.  As  an  initial  step, 
the  ratio  of  v ^/v  is  assumed  to  be  unity.  At  high  wall  temperatures,  the 
exponential  term  will  decay  much  more  slowly  with  y.  For  example,  on  a  wall 
with  a  temperature  of  <jC°F  at  sea  level,  A^t  =  2.25  *  10~\  On  a  wall  with 
a  temperature  of  300Cr  i  at  50, COO  feet,  A  Vu  =  5*46  x  10  or  more  than 
20  times  as  large. 


6.2.2  Viscosity  in  the  outer  region. 

In  the  work  on  incompressible  flows,  the  form  for  eddy  viscosity  in  the 
outer  region  suggested  ty  Clausor  [5]  was  used;  that  is, 
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eo  =  k2P  Ue5*  (6.25) 

where  the  constant  is  taken  to  be  0.0168,  the  value  given  in  Tl8].  The 

same  formula  is  used  for  compressible  flow,  except  that  p  is  a  variable; 
that  is,  p  =  p(x,  y).  Equation  (6.25)  is  modified  by  an  intermittency 
factor  7  that  was  obtained  by  Klebanoff  [19].  It  is  given  by 

7  =  \  [1  -  erf  5  (J  -  .78)]  (6.26) 


where  6  is  the  thickness  of  the  boundary  layer.  This  formula  was  deduced 
from  measurements  of  an  incompressible  flow,  and,  for  want  of  anything  better, 
it  is  also  used  for  compressible  flow.  With  (6.26),  the  eddy-viscosity  formula 
for  the  outer  region  becomes 

eo  =  p  k2  Ue  '*  7  (6.27) 


6.2.3  Definition  of  inner  and  outer  regions. 

The  constraint  used  to  define  the  end  of  the  inner  region  and  the  be¬ 
ginning  of  the  outer  region  is  the  continuity  of  the  eddy  viscosity.  It  can  be 
seen  from  (6.2h)  and  (6.27)  that,  at  a  given  position  along  the  body,  in¬ 

creases  with  y  and  eQ  remains  constant  over  practically  the  whole  boundary 
layer.  Hence,  from  the  wall  outward,  the  expression  for  inner  eddy  viscosity 
applies  until 


c 


i 


=  e 

o 


(6.28) 


or,  in  terms  of  the  distance  from  the  wall,  inner  and  outer  regions  thus  can  be 
defined  as 


c 


i 


0  «  y  5  yc 


(6.29) 


€ 

O 


p  ko 
2  e 


6  7 


yc  ^  y  s  8 


where  yc  is  determined  by  (6.28). 


(6.30) 
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6.2.4  Turbulent  Prandtl  number. 

The  turbulent  Prandtl  number  is  a  measure  of  the  ratio  of  eddy  viscosity 
to  eddy  conductivity,  that  is,  the  ratio  of  the  transport  of  momentum  to  the 
transport  of  heat.  Since  the  flow  in  the  outer  region  shows  some  similarity 
to  a  wake  flow,  one  may  argue  that  a  more  realistic  formulation  of  turbulent 
Prandtl  number  requires  a  separate  expression  in  each  region,  as  in  eddy-vis¬ 
cosity  formulation.  The  fact  that  in  boundary  layers  the  ratio  of  eddy  con¬ 
ductivity  to  eddy  viscosity  is  smaller  than  in  free  turbulence  permits  one  to 
conclude  that  the  lowering  of  this  ratio  is  due  to  the  influence  of  the  wall 
[20].  Consequently,  if  the  ratio  of  eddy  conductivity  to  eddy  viscosity  is 
lowered  by  the  effect  of  the  wall,  it  follows  that  this  ratio  decreases  with 
decreasing  wall  distance  and  increases  with  increasing  wall  distance.  It 
appears  that  at  large  wall  distance  this  ratio  approaches  the  value  2,  that  is, 
the  same  value  observed  for  free  turbulence.  On  the  other  hand,  no  experimental 
results  have  been  obtained  on  the  minimum  value  of  this  ratio  in  the  immediate 
neighborhood  of  the  wall.  For  these  reasons,  as  an  initial  step  the  turbulent 
ftandtl  number  is  assumed  to  be  a  constant  and  equal  to  unity. 

6.3  Transformation  of  Boundary-Layer  Equations 

Before  (6.13),  (6.l4),  and  (6.15)  can  be  solved,  by  a  method  to  be  des¬ 
cribed  later,  it  is  convenient  to  transform  them  to  a  coordinate  system  that 
removes  the  singularity  at  x  =  0  and  stretches  the  coordinate  normal  to  the 
flow  direction,  as  is  usually  done  in  laminar  flow.  First,  the  equations  are 
placed  in  am  almost  two-dimensioneQ.  form  by  the  Probstein  and  Elliot  transfor¬ 
mation  [21]. 


PROBSTEIN -ELLIOT  TRANSFORMATION 


(6.31) 

(6.32) 


where  rQ(x)  is  specified  by  the  body  shape  and  r(x,  y)  is  given  by  (see 
figure  4) 


4 

Figure  4. 'Coordinates  (or  axially  symmetric  body. 


r(x,  y)  =  r0(x)  +  y  cos  a 


(6.1 


Then  from  (6.3 l)  and  (6.3 2), 


d  _  dx  d  +  d 

ox  5x  5?  Sx 
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(6.3 

(6.3 


Define  a  stream  function  f  that  satisfies  the  continuity  equation  (6.13) 
namely, 

di  k 

37  =  r  p  u  (6.36a) 


k/ 

r  (p 


v  + 


p'v* 


6b} 


(6. 
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Now  if  in  the  barred  plane  the  continuity  equation  (6.13)  is  written  as 


(p  u)  +  -^  (p  v  +  p*v* )  =  0 


(6.38) 


then  the  stream  function  ♦  that  satisfies  this  equation  is 


d*  — 

— 1  =  p  u 

dy 


^4  =  -  (p  V  +  P*V* ) 

dx 


(6.39) 


Therefore,  if  the  relations  defined  by  (6. 3M>  (6.35)  (6.37)  and  (6. 39)  are 
used,  (6.36a)  and  (6.36b)  become 


L  /  _r^  dj  A 
W  rk  \  L  *7  /  *7 


p  u  = 


dy  7  dy 


P  v  +  p 


2k 

§  --M-J  ^+^l) 

r  r  v  L  dx  dy  ' 


(6.4oa) 


(6.4ob) 


By  substituting  from  (6.40)  into  (6.l4)  and  using  the  relations  given  by 
(6.34),  (6.35)>  and  (6.39)>  the  following  Probste in -Elliott  transformed 
momentum  equation  is  obtained: 

^V€-=(“^N)  =  -^  +  4[d  +  t)2Vd  + 1  )4(^s )] 


dy  dx  ^  P  dy  ^  dx  dy  dy 


dx  dy 


dy  M  dy 


(6.M) 


where  t  is  the  transverse -curvature  term  defired  as 

t  =  cos  a  (6.42) 

r 

o 

by  using  the  relationship  between  r(x,  y)  and  r  (x)  given  by  (6.33). 
Note  tliat  in  (6.4l)  the  eddy-viscosity  term  has  a  bar.  This  is  necessary 
because  e  is  not  a  scalar  function  (e.g.,  like  p)  and  must  be  transformed 
by  (6.31)  and  (6.32). 


Similarly,  the  energy  equation  (6.15)  may  be  transformed  into  the  form 


I 


!*2h_2I^U4[u  +  t)2N£r(i  + 

iy  Sx  8x  dy  dy  ^  ^  V  p  PrT  '  dy  P 

.  21  A.  ( 1  Al  U  ] 

dy  dy  'p  dy  /y  ^ 


(6.43) 


by  using  (6,4o),  (6.34),  (6.35),  and  (6.39).  Equations  (6.4l)  and  (6.43) 
have  the  desired  two -dimensioned,  form. 

Next,  the  Levy-Lees  transformation  [22]  is  introduced,  in  order  to  put 
(6.4l)  and  (6.43)  into  a  still  more  convenient  form. 


LEVY-LEES  TRANSFORMATION 


p  |i  u  dx 
e  e  e 


(6.44) 


>1/2^ 


(6.1*5) 


A  dimensionless  stream  function  f  is  rexated  to  t  as  follows: 


*  =  (20  '  f(l,  n) 


(6.46) 


With  (6.44),  (6.45),  and  (6.46),  the  partial -derivative  operators  in  the  new 
coordinate  system  ({,  tj)  become 


=  P  ^ 
e  e 


Ue(&  +  St  ^i) 


(6.U7) 


(2fc)l/2  ^ 


(6.48) 


%  -  ^  %  ♦  sf 


(6.49) 


i  -  <2«>1/2  f 


(6.50) 
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where  the  prime  on  f  denotes  the  derivative  with  respect  to  tj.  Introducing 
the  transformations  from  (6.4 7)  through  (6. * iO)  into  (6.4l)  and  (6.4 3)  and 
using  Euler's  equation,  namely, 

dp  +  p  u  du  =0 
e  e  e 


to  replace  the  dp/dx-term  in  (6.4l)  yields  the  transform*:  d  momentum  and  energy 
equations  for  the  compressible  turbulent  boundary  layer. 


MOMENTUM 

or  -|  1  p  -1 

(1 + T)2kc(i  +  ^ )  f  +  ff  +  p  -  (f)2  =25  f  -  f  || 


(6.51) 


(6.52) 


where 


and 


2i 

u  d$ 


C  =  £Ji- 
pePe 


(6.53) 

(6.^4) 


The  prime  on  g  denotes  the  derivative  with  respect  to  and  the  sign  ~ 
indicates  that  the  quantity  is  transformed  by  (6.44)  and  (6.45).  For  example, 
the  transformed  transverse-curvature  term  in  (6. 5l)  and  (6.52)  is  given  by 


(6.55) 


In  (6.51)  and  (6.52)  the  relationship  between  the  quantities  in  the  trans¬ 
formed  plane  and  those  in  the  physical  plane  is  given  by  (6.5 l),  (6.52), 

(6.44),  and  (6.45)  —  what  might  be  called  a  Probstein-Elliot-Levy-Lees  trans¬ 


formation. 
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PROBSTE  IN -ELLIOT -LEVY -LEES  TRANSFORMATION 


d!  =  pelleue(ro/L)2 


(6.56) 


"1/2  L 


(6.57) 


It  can  be  seen  from  (6. 51)  and  (6.52)  that  setting  k  =  0  reduces  them  to  two- 
dimensional  form.  On  the  other  hand,  for  axisymmetric  flow  with  no  transverse- 
curvature  (TVC)  effect,  k  =  1  and  t  =  0,  which  indicates  that  the  ratio  of 
r  to  rQ  is  unity. 


The  boundary  conditions  of  (6.51)  and  (6. 52)  can  now  be  transformed.  From 
(6.1+8)  and  (6.50), 

u  =  u  f'  (6.58) 


and  from  (6.1+0b),  (6.1+7),  and  (6.1+6), 


f(s’ 


e  e  e  r 


since  at  the  wall  u  =  0,  p'v*  =  0,  and  r  =  r^The  length  L  is  taken  oqual  to 
unity  for  two-dimensional  and  for  axisymretric  flow  with  no  TVC-effect.  Then 
in  the  |,  T)-plane  the  boundary  conditions  given  by  (6.19)  become 

1  ^  1  pwv  - 

f({,  0)  =  f  =0  or  f  - r-/o  /  — 7-  - 7—  d|  (mass  transfer) 

w  ”  (se)1/2  y  rk  We 

°  °  (6.55a) 


f’U,  0)  =  0 


(6.59b) 


lim  f'(t,  B)  =  1 

T)  — *  CD 


(6.59c) 


Similarly,  from  the  definition  of  g,  namely,  g  =  H/H  ,  the  boundary  con¬ 


ditions  given  by  (6.20)  become 


lU,°)  =  -ff  or  °) 


(6.60a) 


liin  n)  =  1 

T)  -*  QO 


(6. 60b) 


6.U  Trans  format  ion  of  Eddy -Viscosity  Equations 


The  eddy-viscosity  term  appearing  in  the  momentum  equation  (6. 51)  is  a 
function  of  flow-field  quantities  and  must  be  transformed.  The  turbulent 
Prandtl  number,  on  the  other  hand,  is  assumed  to  be  a  constant  and  therefore 
is  not  transformed. 


Using  (6.3 5)>  (6.Uoa),  (6.^8),  (6.50),  and  the  transformation  given  by  (6.5 7) 
transforms  the  expression  for  the  inner  eddy  viscosity  given  by  (6.2b)  to 

V(p/pe)2kl(21)  /  rkf"  (J~f  J  d’>)  £  '  /  °A  ~k-  ' 

•i - —  J 

(6.6 1) 

Using  (6. 58),  the  relation  given  by  (6.57),  and  the  definition  of  dis¬ 
placement  thickness,  namely, 


00 


d> 


6*  /  (1_fijs_)ay.lEDi^s!/a_(p/p-f.)d, 

J  PeV  “e  °e  y  rk  e 


(6.62) 


transforms  the  expression  for  the  outer  eddy-viscosity  given  by  (0.27)  to 
<„  -  0.00W  -fi-  (201/2l[i  -  erf  5  J  \  di)  -  0.  7s} 

e  n  r  '  n  r 


OO 


0 


(Pc/P  “ 
r 


(6.63) 


6*5  Fluid  Properties 


Fluid  properties  that  appear  in  the  momentum  and  energy  equations  are 
density  (p),  viscosity  (n),  specific  heat  at  constant  pressure  (cp),  and  ther¬ 
mal  conductivity  (k).  The  la4  ,er  appears  in  the  energy  equation  through  the 
laminar  Prandtl  number,  Pr,  defined  as 
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M.  C 


(6.64) 


These  fluid  properties,  which  are  assumed  to  be  functions  of  temperature  only, 
are  given  by  the  following  formulas  [23]  for  air: 


5 


CP(T)  -h  ViTi 

(6.65) 

i=0 

5 

p(T)  «  'L 

(6.66) 

i-0 

Mt)  =  p(cp  +  £  R) 

(6.67) 

where  R  =  1716  lb„-ft/lb  -deg  R  and  the  A's  are  different  for  c  and  n. 

I  m  p 

By  means  of  (6.63),  specific  enthalpy  cam  be  expressed  in  the  form 


h(T) 


6 

L 

i=l 


(6.68) 


This  relation  is  necessary  because,  once  a  solution  of  the  energy  equation  is 
obtained,  the  temperatures  are  required  in  the  calculation  of  the  new  fluid 
properties.  The  coefficients  of  these  polynomials  are  given  in  Table  I. 


In  addition  to  the  relations  given  by  (6.65)  through  (6.68)  it  is  also 
necessary  to  have  the  relation  of  density  to  temperature.  This  can  be  obtained 
from  the  equation  of  state  and  from  the  assumption  that  static  pressure  re¬ 
mains  constant  within  the  boundary  layer. 


P 


Je 

R  T 


(6.69) 


6.6  Overall  Method  of  Solution 


In  previous  studies [7  through  13]  the  method  used  to  solve  the  momentum 
and  the  energy  equations  was  based  on  ideas  originated  by  Hartree  and  Womersley 
[24].  It  consists  of  replacing  the  {-derivatives  by  finite -difference  relations 


-33- 


while  retaining  the  rj-derivatives.  In  this  way,  &  partial  differential  equation 
at  a  given  {-location  is  reduced  to  an  ordinary  differential  equation  that  can 
then  be  solved  by  various  integration  techniques  or  by  finite -difference  tech¬ 
niques,  in  either  nonlinear  or  linear  form.  In  the  previous  studies,  this  ap¬ 
proach  is  used  quite  successfully  in  the  solution  of  the  equations  of  both  in¬ 
compressible  and  compressible  boundary  layers  by  an  integration  technique. 

The  momentum  equation  is  solved  in  its  nonlinear  form. 

When  the  same  method  is  extended  to  the  solution  of  the  equations  of  the 
incompressible  turbulent  boundary  layer,  the  two  chief  disadvantages  of  the 
method  become  quite  pronounced.  First,  computation  time  is  too  long.  This  is 
due  to  the  greater  thickness  of  the  turbulent  boundary  layer  and  to  the  much 
greater  variation  in  the  transformed  velocity  gradient  of  the  turbulent  boundary 
layer.  For  example,  for  a  laminar  boundary  layer  the  transformed  boundary- 
layer  thickness,  t^,  is  about  6,  and  the  transformed  velocity  gradient  at  the 
wall,  f”,  varies  between  0  and  2.  On  the  other  hand,  for  a  turbulent  boundary 
layer  the  transformed  boundary-layer  thickness  may  be  150  or  more,  and  the 
transfrrmed  velocity  gradient  at  the  wall  may  vary  between  0  and  50  or  more. 

The  old  shooting  method,  which  is  based  on  the  cut-and-try  technique,  thus  has 
a  much  bigger  range  to  search  in  meeting  the  boundary  conditions.  Second, 
accuracy  is  reduced  when  small  step  sizes  are  used  in  the  streamwise  direction. 
This  is  a  serious  problem,  because  short  steps  are  essential  near  separation 
or  in  any  region  in  which  changes  in  the  outer  flow  take  place  rapidly. 

To  overcome  both  of  these  disadvantages,  the  momentum  equation  is  linearized 
and  the  integration  method  is  replaced  by  an  implicit  finite -difference  method 
containing  a  variable  grid  in  the  ri-direction.  The  previous  technique  of  re¬ 
placing  the  {-derivatives  by  finite -difference  relations  is  retained.  Solution 
of  the  resulting  algebraic  equations  in  both  momentum  and  energy  equations  is 
obtained  by  the  Choleski  matrix  method  [25]. 


-34- 


6.6.1  Solution  of  the  momentum  equation. 

6.6. 1.1  Finite -difference  representation  of  5 -derivatives.  With  three- 
point-finite -difference  formulas  for  the  {-derivatives  at  £  =  *n  (6. 5l) 
becomes 


(6.70) 


=  25, 


n 


f' 

n 


A  f  •  +  A  f '  .+  A,f*  _ 
In  2  n-1  3  n-2 


*  V„-l+  Vn-2 


where  the  coefficients  A^,  A^,  and  A^  are  as  foll(  ws : 


For  three  points : 


Al  '  ~\.l>  +  (S„  ~\.2> 

.  ^n  ~  Sn-2 _ 

2  “  *n-l^*n-l  ~  5n-2^ 

tn  “  in.x 

Aj  =  (Sn  -  Sn-2,({n-l  "  W 


(6.71a) 

(6.71b) 


(6.71c) 


At  5  =  5^  the  quantities  A^,  A^,  and  A^  are  known,  and  the  quantities 
having  the  subscripts  n-1  and  n-2  are  known  functions  of  tj  from  solutions 
obtained  at  the  two  previous  stations.  Thus,  at  5  =  5n>  (6. 50 )  is  an  ordinary 
differential  equation  in  tj.  There  is  no  problem  of  starting  the  solution, 
because  the  terms  with  {-derivatives  disappear  since  {  -  0.  At  the  next  sta¬ 
tion,  5  ,  the  three -point  formulas  are  replaced  by  two-point  formulas;  at  all 
stations  farther  downstream  the  three -point  formulas  are  used. 


For  two  points : 


A  = 

1  "  i  -  I 


A«  =  - 


n-1 


5  -  $  , 

n  n-1 


A^  =  0 


(6.72) 


6. 6. 1.2  Linearization  cf  momentum  equation.  As  is  stands,  (6.  70)  is  an 
ordinary  nonlinear  differential  equation  of  third  order,  and  with  boundary 
conditions  given  by  (6.59)  it  is  difficult  to  solve.  In  [7  through  13]  this 
equation  is  solved  :juite  satisfactorily  for  laminar  and  incompressible  turbu¬ 
lent  boundary-layer  flows  in  its  nonlinear  form;  here  it  is  solved  in  a  linear 


i 
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form  in  order  to  speed  up  the  computations  and  to  be  able  to  take  smaller  step 
sizes  in  the  streamwise  direction. 

Two  different  linearization  techniques  were  used  to  linearize  (6.70).  The 
first  one  is  based  on  the  quasi -linearization  technique  [26].  If  (6.70)  is 
written  in  the  form 

g(f,  f',  f",  f"*,  TJ,  0  =  o  (6.73) 

and  expanded  in  a  Taylor  series  around  a  known  solution  designated  by  subscript 


the  following  quasi-linearized  equation  is  obtained: 

E1qp,M  +  E29"  +  E5<p»  +  E^qp  =  E5  (6.75) 

where  qp's  are  the  perturbation  quantities  defined  as 


<p  =  f  -  f  ,  cp*  =  f*  —  f',  etc. 
o  o 

Note  that  in  previous  studies,  <p*s  are  used  to  denote  translated  stream 
function  f's  and  should  not  be  confused  with  the  perturbation  terms  that  are 
denoted  by  <p*s  in  this  study.  The  latter  are  used  for  what  are  commonly  called 
e's,  because  e's  have  already  been  used  to  denote  eddy-viscosity  terms. 


The  coefficients  Ej  are  given  by 

E  =  (1  +  T)2k  C  (l  +  -  ) 

1  ov  4  'o 

Pk  ^  I 

=  [(1  *  t)  Co(i  .  -  )o]  .  f.  2S(AlV  A2fn.1+  A5fn.2) 

EJ  -  -  2(0  @  ♦  t(2Alf*  .  A2f-.1  .  A5f'.2)] 


(6.76a) 
(6. 76b) 


(6.76c) 
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E4  =  f"  (1  +  2|  kj  (6. 76d) 

e5  .  -  d  ♦  t)2k  cod  *  *  +  t)a  cod  +  £  )o]'  r-  -  y» 

-  @  t(oe/p)0  -  (f-)2]  -  2irr;<v;  +  Vi-1  *  Vi-2> 

-fo'Vo+ Vn-l*  Vn-2»  ««•«■) 

When  (6.75)  is  solved  by  the  finite -difference  method  to  be  described  later, 
successive  iterations  on  f”  cause  oscillations  in  f”  and  convergence  is 
not  obtained?  On  the  other  hand,  (6.75)  shows  no  oscillations  in  f£  for 
laminar  flow.  For  the  latter  case  the  convergence  is  fast  (quadratic),  and  the 
results  are  very  satisfactory.  However,  unlike  the  old  shooting  method  used  in 
previous  studies,  this  method  gives  results  that  indicate  that  mu^h  finer 
spacing  in  the  strearawise  direction  is  necessary  to  obtain  results  comparable 
to  those  of  the  old  shooting  method.  That  is  to  be  expected,  since  with  finer 
spacing  the  difference  between  the  calculated  and  assumed  solutions  (f  and  its 
derivatives)  decreases. 

The  second  technique  used  to  linearize  (6. 70 )  is  in  principle  similar  to 
the  one  used  in  the  energy  equation  in  previous  studies  (see  [11], for  example). 
In  this  case,  certain  terms  that  make  the  equation  nonlinear  are  assumed  to  be 
known  from  previous  iterations,  that  is. 


Since  the  previous  studies  have  indicated  that  the  greatest  error  usually 
appears  in  f^  ,  again  successive  values  of  f£  are  used  for  the  convergence 
criterion;  that  is, 


where 


|j»n  _  ^ti 

V  V 

Vl  Q 

7  is  a  prescribed  value  and  Q 


<  7 

denotes  the  iteration  number. 


[(1  ♦  t)2kCo(l  ♦  J  )0f]’  *  f„f"  *  P[(ae/0)o  -  fof'J  ■ 

«[  w ♦  v;-x  *  v;.2)  -  v  (Aif  +  v.**  v-p] 

The  subscript  o  denotes  that  the  function  is  obtained  from  a  previous  itera¬ 
tion.  For  simplicity,  the  subscript  n  is  dropped.  Equation  (6.77)  can  be 
expressed  in  the  same  form  as  (6. 75),  where  again  the  <p's  are  the  perturbation 
quantities  defined  as  before.  However,  the  coefficients  Ej  now  are  given  by 

E.  =  (l  +  T)2k  C  (l  +  -  ) 

1  o  n  7o 


E  =  (l  +  ?)2k  C  (1  +  -  )  ♦  f 

2  'ono  o 


E,  =  -  Pf*  -  2{f'A.  =  -  f*  (3  +  2{  A  ) 
3  o  o  1  o  1' 


(6.78) 


E4  =  2t  f*«Al 

E5  =  -  3(pe/p)o  ♦  2«f-(Vi.1  ♦  A3f*.2)-  f”(A2fn_1+  A5fn.2)] 

-  [E  f  •  *  *  +  E„f"  +  E,f»  +  E,  f  ] 

1  o  2  o  3o  4oJ 

A^ain,  when  (6.75)  with  the  coefficients  defined  by  (6.78)  is  solved  by 
the  same  finite-difference  method  as  before,  oscillations  in  f"  continue 
and  convergence  is  not  obtained.  On  the  other  hand,  the  results  for  laminar 
flow  indicate  that  even  though  the  convergence  rate  of  this  linearization  method 
is  not  as  rapid  as  the  other  one,  this  linearization  methods  shows  much  less 
sensitivity  to  streamvise  spacing  and  the  results  are  as  accurate  as  the  other 
one.  Furthermore,  when  a  3-point  mean  of  the  eddy-viscosity  expression  appear¬ 
ing  in  (6.78)  is  taken  in  the  n-direction  and  a  2-point  mean  of  the  coefficients 
given  by  (6.78)  is  taken  in  the  {-direction,  oscillations  in  f£  stop  and  ex¬ 
cellent  convergence  on  f”  is  obtained. 

6, 6. 1,3  Variable-grid  spacing  in  the  redirection.  From  a  computational 
aspect,  a  turbulent  boundary  layer  presents  a  much  more  difficult  problem  of 
calculation  than  a  laminar  boundary  layer.  Consider,  for  example,  an 
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incompressible  turbulent  flow.  The  skin-friction  is  appreciably  greater  than 
it  is  for  laminar  flow.  This  means  that  du/dy  is  greater.  To  maintain  com¬ 
puting  accuracy  when  du/dy  is  large,  short  steps  in  y  must  be  taken;  v/hen 
it  is  small,  longer  steps  can  be  taken.  Therefore,  near  the  wall  the  steps 
in  a  turbulent  boundary  layer  must  be  shorter  than  they  are  in  a  laminar  layer 
under  similar  conditions.  In  the  outer  region,  the  expression  for  the  eddy 
viscosity  can  be  written  as  (see  6.27) 

=  0.0168  (p/p  )  R  #  (6.79) 

e  6 

#  . 

where  R  ,  is  defined  as  p  u  6  /p  and  y  is  assumed  to  be  unity.  Since 

R  #  may  easily  exceed  10,000,  e  /p  may  easily  exceed  100.  This  part  of  tiie 

^  o  e 

boundary  layer  then  behaves  like  that  of  a  highly  viscous  flow  (oil,  for  example), 
and  the  thickness  of  this  part  becomes  very  great.  On  the  whole,  the  turbulent 
boundary  layer  is  something  like  a  cold  oil  flowing  past  a  heated  wall.  Thus 
it  can  be  seen  that  the  turbulent  boundary  layer  is  characterized  by  great  thick¬ 
ness;  but,  in  spite  of  that,  the  steps  near  the  wall  must  be  shorter  than  for 
laminar  flow.  As  a  result,  if  steps  of  constant  q-spacing  are  used,  far  more 
are  required  than  are  necessary  to  solve  a  laminar  layer,  nor  will  ordinary 
b^ing  improve  the  situation. 

A  possible  solution  of  this  problem  is  to  devise  a  simple  variable -step 
system  that  has  short  steps  near  the  wall,  which  lengthen  with  distance  from 
the  wall.  A  promising  idea,  illustrated  in  figure  5>is  a  grid  whose  spacing 
is  such  that  the  ratio  of  lengths  of  any  two  adjacent  intervals  is  a  constant; 
that  is,  h^  =  Kh^  The  distance  to  the  i-th  grid  line  is  given  by  the 
following  formula: 

\  =  hiy~K~'-~'l)  1  =  °»  1»  2*  3,  ...  N  (6.80) 

There  are  two  parameters:  h^,  the  length  of  the  first  step,  and  K,  the  ratio 
of  two  successive  steps.  A  number  of  useful  relations  for  this  system  can  be 
derived,  but  the  following  is  of  particular  interest: 


(6.81) 
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Figure  5.- Finite-difference  variable-grid  system  in  the  q-direction. 

Suppose  the  boundary  layer  has  a  thickness  =  100.  Assume  that  an  error 

analysis  shows  that  h-  must  be  0.1,  but  that  it  also  shows  that  h  (=  h  ) 

X  i '  m' 

could  be  2.  Wl.at  is  K?  It  is 


100  -  0.1 
100-2 


^  =  1.019^ 


and  157  steps  are  required  for  the  traverse.  At  a  constant  step  length  of  0.1, 
which  is  based  on  the  wall  accuracy  requirements,  1000  steps  would  be  required. 

In  the  usual  problem,  the  value  of  K  lies  between  1.01  and  1.02.  Figure  5 
accurately  represents  the  spacing  for  K  =  1.07,  a  grossly  large  value.  The  spa¬ 
cing  in  the  {-direction  is  arbitrary. 

6.6.1.J*  Method  of  solution.  An  implicit  finite-difference  technique  is 
used  to  solve  (6.75).  Figure  6  shows  the  finite-difference  molecule  selected. 

The  variable-grid  differentiation  formulas  for  this  molecule,  as  well  as  for  a 
3-point  molecule  that  is  used  in  the  energy  equation,  are  de -ived  for  the  Ta grange 
interpolation  formula  and  are  given  in  Appendix  C  (see  also  [27]).  By  using 
the  first-,  second-,  and  third-derivative  formulas  at  point  (n,  i),  (6.75)  can 
be  written  as  follows : 


all-1 


-4i- 


Ri  =  B5 


E  E  "| 

6(ax+  K  -  K2)  +  2  ^  hi-2K(ai“  ^1“  K2)  “  ^  hi-2K5ai  J  (6.8?f) 


The  eddy-viscosity  expressions  appearing  in  the  coefficients  E^  and  E^  at 
(n,  i)  are  taken  as  the  J -point  means  of  the  eddy-viscosity  expressions  in  the 
rj-direction;  that  is. 


>  =  ± 


av'i 


i-1 


+  e.  +  c 


i+1- 


(6.84) 


This  averaging  process  was  introduced  in  order  to  stop  the  oscillations  in  f^ 
and  cause  the  iterations  to  converge. 


For  each  point  (n,  i)  at  station  i,  an  algebraic  equation  of  the  form  of 
(6.82)  is  written,  yielding  N  —  3  algebraic  equations  with  N  unknowns.  The 
other  three  equations  are  obtained  by  considering  the  following  finite -differ¬ 
ence  molecule  at  (n,  N-l)  (see  figure  7)  and  by  the  boundary  conditions  given 
by  (6.39). 


5 

Figure  7  -  Finite-difference  molecule  for  the  momentum  equation  at  (n,  N-1). 


By  using  the  variable-grid  differentiation  formulas  at  pr  int  (n,  N-l),  (6. 73) 
can  be  written  as  follows: 


%  +  FN-1%-1  +  GN-l%-2  +  In-l%-3  +  ^-l^-1*  =  NN-1 


(6.85) 


-42- 


where 


B, 


N-1  Vj 


E 

6(a2+  Ka.+  K2-  K3)+  2g£  h^K^a^  Kag-  K*a2+  K^-  k\±-  K4) 

E  E  "I 

+  E?  hW-4K  ^aia2”  K  aia2”  k2&2~  K\)  +  B,E  hu  ,  J  (6*86a) 


N-l 


B, 


RN- 


6(a2+  K  ar  K3)+  2  ^  h^K^ 


B4ElhN-4 

-  K2^-  K^) 

E 

-  if  ^-4^  aia2^ 


(6.86b) 


B„ 


hi-l  ~ 


H-l 


E  E 

fi( V  K2-  k3)+  2  e7  hM-4K2(a2-  K  a2  -  k5>-  *4j,(k5*2> 


(6.86c) 


B  jj  £ 

3.-1  =  5^  i6K<al  +  K  -  K2)*  2  ^  K  V  «?)  -  ^  hN.lA 


(6.86d) 


N, 


N"1  ElRN-lhN-4 


(6.86e) 


*H-1  ■  B5 


E  E 

6<V  Ka1+  K2)  +  2  ^  ^(Ka^  ^2+  k\)+  g  h^K^a, 


(6.86f) 


By  using  the  first -derivative  variable -grid  differentiation  formula  at  point 
(n,  0),  the  boundary  condition  given  by  (6.59b)  can  be  written  as  follows: 


where 


%  +  Fl^  +  Gl*2  +  =  N1 

b4 

F1  =  ^  ala5 

B3 

G1  =  a^a3 


(6.87) 

(6.88a) 


(6.88b) 


(6.88c) 
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Fl  = 


B2 

ala2a3 
B 


Nl  = 


^aia2+  aia3+  a2a3+  aia2a3^fw  “  [Il^fo^l  + 
Gl<fo>  +  F1  (fo>3+  (fo>J 


=  B^a,a 
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(6.88d) 

(6.88e) 


The  term  f  is  given  by  (6.59a).  It  is  equal  to  zero  when  there  is  no  mass 
transfer. 


Similarly,  by  using  the  first-derivative  variable -grid  differentiation 
formula  at  point  (n,  N),  the  boundary  condition  given  by  (6. 59c)  can  be  written 
as  follows: 

%  +  FN%-1+  GN%-2+  HNC%-3+  IN<%-1*  =  NN  (b.89) 


where 


b4  3 

FN  Rjj  K  aia2a3 


(6. 90a) 


B3  4 
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(6.90b) 


^  Rjt  K^ala3 


(6.90c) 
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Rjj  K  aia2 


(6.90d) 


Nn  w>  r  AAVl*  GN^fo^N-2+  Wn-3+  Vfo^N-4+  AV  (6*90e) 


Rjj  =  B5(K5a1a2a5+  K^ay-  K^a^  K6a1a2) 


(6.90f) 


-1+1+ - 


At  station  n,  (6.82),  (6.85),  (6.87),  and  (6. 89) form  a  system  of  N  linear 
simultaneous  algebraic  equations  with  N  unknowns  that  can  be  represented  by 
a  matrix  equation  of  the  form  A  <p  =  N. 


Fl+  1 

•  •  • 

•  •  •  • 

^-3  ^-3  GN-3  FN-3  -1 

”n-2  FN-2  °N-2  FN-2  1 

^-1  ^-1  °N-1  FN-1  1 

^  ^  gn  fn  1 


(6.91) 

The  Choleski  method  [25]  is  used  to  solve  (6.9l)»  The  matrix  A  is  expressed 
as  the  product  of  a  lower  triangular  matrix  L  and  an  upper  triangular  matrix 
U;  that  is 


-1*5- 


so  that  (6.91)  becomes 


L  U  <p  =  N 

(6.93) 

If  (6.93) 

is  written  in  the  form 

L  Y  =  N 

(6.94) 

where 

U  q>  -  Y 

(6.95) 

then  the  method  of  solution  is  apparent.  The  two  triangular  matrices  L  and 
U  are  found  from  (6.92),  the  column  matrix  Y  from  (6.9I*),  and  the  column 
matrix  q?  from  (6.95)* 


6.6.2  Solution  of  the  energy  equation. 

6. 6.2.1  Finite-difference  representation  of  5-derivatives.  The  method 
of  solution  of  the  ener gy  equation  is  similar  to  that  of  the  momentum  equation. 
Arain  the  S-derivatives  are  replaced  by  finite  differences  that  are  defined 
by  (6.71).  Thus  (6.52)  is  written  as 


lp;W  +  (VA  =  (p2„ 


(6.90) 


where  t'  "'efficients  P,  are  given  Ly 

J 


rT 


(6.97a) 


9  "  f‘  -  25(V  *  Vn-i  *  V„-2 


(6.97b) 


Pv  =  -  25  f‘  Ax 


(6.97c) 


h  =  -6  ^(A^  +  A7;gn_J  -1(1+  t)2k  C  (l  -  Yp)f*Cn 


(6. 97d) 


(6.99d) 
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Again,  for  simplicity,  the  subscript  n  is  dropped.  In  solving  (6.98),  values 
of  f  and  its  derivatives  determined  from  the  previous  solution  of  the  momen¬ 
tum  equation  and  fluid  properties  determined  from  the  previous  solution  of  the 
energy  equation  are  used.  Details  of  the  iterative  procedures  for  solving  both 
momentum  and  energy  equations  are  given  in  Section  6.6.5.  With  these  proce¬ 
dures,  (6.98)  is  linear  and  g  is  tiie  only  unknown. 

For  each  point  (n,  i),  at  station  i  an  algebraic  t?quation  of  the  form 
of  (6.98)  is  written,  yielding  N-l  algebraic  equations  with  N-l  unknowns. 

Tie  function  g  at  N  is  known  from  (6.60b),and  g  at  the  wall  is  known  from  one 
fern*  of(6.u0a),  depending  on  whether  tlie  temperatuic  or  the  heat  transfer  is 
specified  at  the  wall. 

Case  l;  g^  is  known. 

For  i  =  1,  (6.98)  becones 

g2  +  Slgl  +  Vw  =  Z1  (6.100a) 

and,  since  gw  is  known,  it  can  be  written  as 

g2  +  Slgl  =  V  Wl®w  "  Z1  (6.100b) 

Case  2:  g*  is  known. 

For  this  case  cw  is  unknown,  but  since  is  known,  gw  can  be  ex¬ 

pressed  in  terms  of  g^  and  g^'s  by  using  (6.100a)  and  the  first-derivative 
variable-grid  differentiation  formula  at  point  (n,  0). 

I 

gj  =  “  h1  B1(2  +  K)rv  +  B2(l  +  K)gx  +  B3g2  (6.101) 

Thus,  by  eliminating  gw  from  (6.100a)  and  (6.101)  and  rearranging, the  following 
equation  is  obtained: 


The  system  of  equations  given  by  (6.98)  is  also  represented  by  a  matrix 
equation  similar  to  (6.91),  and 


(6.103) 


is  solved  by  the  Choleski  method  described  in  Section  6.6.1.^.  In  the  first 

row  of  (6.103),  the  asterisks  on  the  coefficients  and  indicate  that 

they  are  determined  by  the  boundary  conditions  at  tlie  wall.  When  n  is 
*  *  " 
specified,  S..  and  Z  are  defined  by  the  coefficients  of  (6.100b),  and  when 

g^  is  specified,  and  Z^  are  defined  by  the  coefficients  of  (b.102). 


6.6.3.  Outline  of  procedure  for  solving  momentum  and  energy  equations  simul- 


The  general  method  of  solving  the  momentum  and  energy  equations  at  station 
n  is  described  in  Sections  6. 6. 1.4  and  6.6.2.?,  respectively.  The  procedure 
for  solving  both  of  these  equations  at  station  n  is  discussed  in  this  section. 
Consider  the  case  when  the  program  is  solving  the  momentum  and  energy  equations 
at  station  n.  Values  of  f  and  g  and  their  derivatives  at  all  previous 
stations  will  be  known.  In  addition,  the  fluid  properties  are  known  up  to 
station  n.  Calculations  start  from  the  momentum  equation.  Before  the  momen¬ 
tum  equation  can  be  solved,  it  is  necessary  to  establish  the  inner  and  outer 
regions.  Since  the  eddy-viscosity  expressions  contain  terms  like  f"  and  6  , 
these  two  regions  are  not  known  until  a  solution  of  the  momentum  equation  is 
generated.  Thus  an  iteration  process  is  necessary.  For  the  first  iteration, 
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i. 


&  and  f"  art  v  btained  from  the  final  solution  at  station  n-1,  and  inner 
and  outer  regions  are  established  by  (6.28).  With  this  information,  the  system 
of  equations  represented  by  (6.91)  is  solved  by  using  the  fluid  properties  from 
the  n-1  station  and  f  and  its  derivatives  from  the  n-2,  n-1  stations.  Once 
(6.91)  is  solved,  new  values  of  f  are  obtained  from  qp*s  and  f's.  The  first 
derivative  of  f  is  subsequently  calculated  from  the  values  of  f,  the  second 
derivatives  from  the  values  of  f',  and  the  third  derivatives  from  the  differen¬ 
tial  equation,  (6.75)*  A  solution  at  station  n  for  one  iteration  is  then 
defined  as  the  averages  of  the  calculated  values  of  f  and  its  derivatives 
and  previous  values  of  f  and  its  derivatives;  that  is,  the  values  of  f  and 
its  derivatives  are  obtained  from 


f  I  I  |‘fc>Q+  (foVl  >  f’  ‘I  +  (f, 


I  ) 

o'Q+lJ’ 


etc. , 


where  Q  is  an  iteration  number.  This  is  necessary  to  stop  the  oscillations 

in  f". 

w 


After  a  solution  of  the  momentum  equation  is  obtained  in  this  manner,  the 
values  of  f  *s  and  their  derivatives  from  this  solution  are  used  to  solve  the 
energy  equation.  Again,  the  fluid  properties  from  station  n-1  are  used. 

On^e  (6.103)  is  solved,  fluid  properties  are  obtained  for  that  particular  solu¬ 
tion  by  converting  the  total  enthalpy  to  static  temperature,  and  inner  and 
outer  regions  for  the  next  solution  of  the  momentum  equation  are  established. 
TJie  conversion  of  enthalpies  to  temperatures  is  necessary  because  the  fluid 
properties  given  by  (6.65),  (6.66),  (6.67),  and  ( 6.6 9)  are  expressed  as  func¬ 
tions  of  temperature.  The  procedure  of  calculating  temperatures  from  the 
enthalpies  is  as  follows:  From  the  steady-state  energy  equation,  the  total 
enthalpy  H  is 


H  I  h  +  |  u2 


(6.104) 


By  using  the  definition  of  g  and  the  transformation  given  by  (6.58),  equa¬ 
tion  (6.104)  can  be  written  as 


g  H  -  h  uf(f')? 


(6.105) 


I 
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Since  g  is  known  from  the  solution  of  the  energy  equation,  f*  is  known  from 

the  solution  of  the  momentum  equation,  and  H  and  u  are  given,  the  right- 

c  c 

hand  side  of  (6.105)  is  known.  Equation  (6.105)  is  effectively  a  sixth-degree 
polynomial  given  by  (6.68).  Rewriting  (6.105)  and  denoting  it  by  F  gives 

F  =  h-  gHe  +  |  ue2(f’)2  =  F(T)  (6.106) 


The  Newton -Raph son  method  is  used  to  solve  (6.106).  An  initial  temperature, 
T^,  is  obtained  from  the  perfect- gas  relationship;  that  is. 


T 


1 


!u;Vr 


where  c 

P 

calcilnJ  1 


^  O 

for  air  is  6.055  x  1(T  ft^/sec-deg 
by  the  formula 


F(T  ) 

t2  =  Ti  ”  FTt^T 


R. 


(6.107) 

A  new  temperature  then  is 

(6.108) 


where  F*  is  the  derivative  of  F  with  respect  to  T.  Once  is  obtained, 

the  sane  procedure  can  be  repeated  and  new  values  of  T  can  be  obtained  until 
the  difference  between  the  consecutive  values  of  T  is  less  than  y,  where  y 
is  a  small  number. 


Once  the  temperature  at  each  point  within  the  boundary  layer  is  known,  the 

fluid  properties  are  calculated  from  (6.65),  (6.66),  (6.67)  and  (6. 69). With  these 

* 

newly  calculated  fluid  properties  and  with  the  values  of  8  and  f  from  the 
previous  solution  of  momentum  equation,  inner  and  outer  regions  are  established 
and  the  momentum  equation  and  the  energy  equation  are  solved  in  succession. 

An  iteration  procedure  based  on  the  convergence  of  f^  is  followed.  The  itera¬ 
tion  process  continues  until  either 


7 


or 


Q  =  Q, 


max 
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where  7  is  a  prescribed  value  ar.d  Q  denotes  the  iteration  number.  A 
typical  value  of  Q,  is  6.  Once  either  of  these  conditions  is  satisfied,  the 
procram  proceeds  to  the  next  station.  Figure  9  shows  the  flow  diagram  at 
station  n. 


Repeat  this  cycle  until 


(V  -  (f;) 

Q+l  Q 


<  7 


or 


Q  "  Snax 


Then  proceed  to  station  (n  +  l). 


Figure  9.- Flow-diagram  for  solving  the  boundary-layer  equations  at  station  n, 


As  the  calculations  proceed  downstream,  the  boundary-layer  thickness  in¬ 
creases.  Since  at  station  n  the  initial  values  of  f  and  its  derivatives 
are  obtained  from  station  n-1,  it  is  necessary  to  make  an  assumption  for  these 
values  for  T]  >  'On-  where  <l>n.  ^  is  the  transformed  boundary-layer 
thickness  at  station  n-1.  For  this  reason,  at  station  n  the  values  of  f 
and  its  derivatives  are  obtained  from  station  n-1  up  to  =  (T\#)n_p*  ^ or 

i  ><!.>».  f  is  obtained  from  f  =  f'  is  assumed 

to  be  unity,  and  f"  and  f”'  are  assumed  to  be  zero  (see  figure  10).  The 
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latter  assumption  is  permissible,  since  f"  and  f'"  approach  zero  as 

n  -*  V 


Figure  10.-Diagram  showing  the  method  of  generating  the  initial  coefficients  of  the  momentum  equation 
at  station  n. 

6.6.4  Starting  the  solution  at  the  leading  edge  or  at  the  stagnation  point, 

Ljl.2- 

At  5  =  0,  the  5-dependent  terms  disappear  in  both  the  momentum  and  energy 
equations.  Hence  the  coefficients  of  (6.75)  given  by  (6.JB)  and  the  coefficients 
of  (6.96)  given  by  (6.97)  do  not  contain  5-dependent  terms.  In  addition,  the 
flow  at  this  zero-station  is  assumed  to  be  laminar,  because  the  eddy-viscosity 
and  turbulent -Prandtl-number  terms  are  zero. 

The  procedure  of  solution  requires  that  the  momentum  equation  be  solved 
first,  but  to  do  so  requires  values  of  the  fluid  properties  as  well  as  the 
initial  values  of  f  and  its  derivatives  denoted  by  the  subscript  o.  The 
latter  values  are  obtained  by  the  old  shooting  method  of  previous  studies 
[7  through  13 J •  The  fluid  properties  are  obtained  by  assuming  a  linear  tempera¬ 
ture  profile  of  the  form 

T  =  a  +  br) 

T  =  T  +  (T  -  T  )-5- 
w  e  w  n 


which  becomes 


(6.109) 


tj  assuming  T(o)  =  Ty,  T(t^)  =  Tg. 


After  the  first  solution  of  the  momentum  equation  is  obtained  by  (6.91), 
new  values  of  f  and  its  derivatives,  together  with  the  fluid  properties  ob¬ 
tained  by  a  linear  temperature  variation,  are  used  to  solve  the  energy  equa¬ 
tion.  Again,  once  (6.10?)  is  solved,  fluid  properties  are  obtained  for  that 
particular  solution  by  converting  the  total  enthalpy  to  static  temperature. 
These  fluid  properties  are  used  for  the  next  solution  of  the  momentum  equation. 
Figure  11  shows  the  flow  diagram  at  5  =  0. 


Repeat  this  cycle  until 


or  4  = 

Then  proceed  to  station  1. 


Figure  1 1 Flow  diagram  for  solving  the  boundary-layer  equations  at  §  =  0. 
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6.7  Boundary-Layer  Farameters 


Once  the  profiles  of  f  and  g  and  their  derivatives  are  determined  at 
any  {-station,  the  boundary-layer  parameters,  such  as  displacement  thickness, 
momentum  thickness,  local  shear-stress  coefficient,  local  skin-friction  co¬ 
efficient,  local  heat-transfer  rate,  Stanton  number,  as  well  as  the  law-of-the 
wall  coordinates,  velocity-defect  coordinates,  and  other  parameters  of  interest, 
can  be  calculated.  Some  of  th*.se  parameters  are  given  by  the  following  equa¬ 
tions. 


Displacement  thickness: 


or,  in  terms  of  {,  q-coordinates, 


# 

6 


u 


Momentum  thickness: 


or,  in  terms  of  5,  rj-coordinates, 


oo 


6  = 


—  J  f'  (l  -  f')  drj 


k  p  u 
r  e  e 
o  0 


Shear  stress  at  the  wall: 


or,  in  terms  of  {,  rj-coordinates 


,W=  I(25)1/2 


p  p  rk  u2  f" 
w  w  o  e  w 


(6.110a) 


(6.110b) 


(6.111a) 


(6.111b) 


(6.112a) 


(6.112b) 
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Local  skin -friction  coefficient: 


c„  = 


f  12 

o  P  ^ 
2  e  e 


(6.115a) 


or,  in  terms  of  |,  Tj-coordinates , 


K. 

i  H  P,., 

:  =  (2/0lA  -r~  —  n  f" 

f  '  L  p  w  w 


(6.115b) 


Local  shear-stress  coefficient  for  laminar  flow,  not  necessarily  at 

tL 

c-tT  1  ? 

L  2  PeUe 


the  wall 

(6.1l4a) 


or,  in  terms  of  £,  ^-coordinates, 


c  =  (2/|) 

tL 


1/2  _L 1 


(6.1l4b) 


Local  shear-stress  coefficient  for  turbulent  flow: 


p  u’v1 

1  2 
o  P  u 
2  e  e 


(6.115a) 


or,  in  terms  of  |,  ^-coordinates, 


^  .  (a/o^-p  S- .  f 


(6.115b) 


Total  shear-stress  coefficient: 


total 


■  +  o  =  (2/6)1/2-S- -£-{(*+  c)f" 

L  T  Me 


(6.116) 


Equation  (6.116)  reduces  to  (6.115b)  at  the  wall,  since  €  =  0  at  y  =  0. 


Heat-transfer  rate  at  the  wall: 


-  s,  *  r-  <»>/*>. 

P,. 


(6.117a) 
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56. 


or,  in  term*  of  ft,  ^-coordinates, 

.Jo  Wo”.  . 

1  Prist??*  “» 


(6.117b) 


Stanton  number: 


St 


p  u  (H  -  H  ) 
e  e  e  w' 


or,  in  terms  of  (,  ^coordinates, 
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St 


ro  pw  % 

L  p.  Prv(2t)1/a(l  -  «.) 


The  lmw-of -the -wall  coordinates: 


“T?  -  (2/cf)1/2f» 


u 

* 

_EL_  = 
v 


(i  01/2l 


Pe  V 


k  (Pe/p)dT 


(6.1l8a) 


(6.118b) 


(6.119a) 

(6.119b) 


The  velocity-defect  coordinate: 


u  -  u 
e 


u 


(2/cf)l/2(l  -  f») 


(6.120) 


Once  &  ,  ©,  and  c^,  are  calculated,  some  of  the  other  boundary-layer  para¬ 
meters  such  as  shape  factor  H,  defect-shape  factor  G,  and  defect-displacement 
thickness  A  can  be  calculated  from  the  following  relations: 


Shape  factor: 


Defect-shape  factor: 


H  =  6  /e 


(6.121) 


00  /  u  —  u  \  2  /  f?  /  u  -  u  \  o\2 


0 


(6.122) 
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Defect-displacement  thickness: 


A 


dy 


(6.123) 
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7.0  COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL  RESULTS 

The  ultimate  test  of  any  numerical  calculation  technique  is  a  comparison 
of  calculations  with  exact  solutions  and  with  experiments.  However,  when  no 
exact  solutions  are  possible,  as,  for  example,  in  turbulent  flows,  all  that 
can  be  done  is  to  compare  the  calculations  with  experiments.  For  this  reason, 
a  variety  of  both  incompressible  and  compressible  turbulent  flows  are  calcu¬ 
lated  by  the  present  method,  and  comparisons  with  experiments  are  made.  The 
results  presented  in  this  report  do  not  represent  a  finished  development,  but 
are  only  what  has  been  obtained  so  far  by  one  particular  formulation  of  eddy 
viscosity  and  turbulent  Prandtl  number. 

In  the  present  method,  the  solution  begins  at  the  leading  edge  or  at  the 
stagnation  point,  where  5  ■  0,  and  proceeds  downstream.  At  station  5=0, 
the  flow  is  laminar,  and  it  becomes  turbulent  at  any  specified  station  where 
5  >  0.  In  principle,  a  calculation  can  be  started  in  the  middle  of  a  flow 
field,  provided  that  information  is  available  to  be  used  as  an  input  from  up¬ 
stream  positions  in  the  flow.  In  some  of  the  experiments,  a  specific  portion 
of  the  flow  was  instrumented.  This  is  referred  to  as  the  test  section. 

Details  of  the  flow  field  upstream  of  the  test  section  were  not  reported,  and 
hence  it  ./as  necessary  to  add,  in  effect,  an  initial  length  (usually  a  flat 
plate)  to  the  test  section.  This  length  was  selected  to  match  as  well  as 
possible  the  calculated  boundary-layer  parameters,  such  as  0  and  H,  with 
the  experimental  values  at  the  beginning  of  the  test  section.  For  each  test 
case,  details  of  the  matching  procedure  are  discussed  separately. 

The  computer  program  of  the  present  method  has  a  capacity  of  500  points 
in  the  q-direction.  There  is  no  restriction  of  the  number  of  stations  in  the 
5-direction.  For  a  turbulent  boundary  layer,  the  transformed  boundary-layer 
thickness  may  be  130  or  more,  and  the  transformed  velocity  gradient  at  the  wall 
may  vary  between  0  and  30  or  more.  Since  a  variable -grid  system  is  being  used 
in  the  rj-direction  and  smaller  spacing  is  necessary  in  the  inner  region,  which 
occupies  approximately  10  to  15  percent  of  the  whole  boundary  layer,  it  is  con¬ 
venient  to  have  a  table  of  variable-grid  parameters  such  as  Table  II.  For  the 
cases  that  have  been  calculated  so  far,  an  initial  rj-spacing  of  h^  =  0.01 
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TAJLE  II 


VALUES  OF  VARIABLE -GRID  PARAMETERS  FOR  DIFFERENT 
VALUES  (N  *  500  POINTS) 


hi  " 

0.001 

it 

0.01 

% 

K 

hedge 

K 

hedge 

30 

1.01196 

0.355 

1.00595 

0.187 

65 

1.01384 

0.888 

1.00807 

0.530 

100 

1.01487 

1.466 

1.00920 

0.922 

130 

1.01553 

2.019 

1.00991 

1.305 

150 

1.01583 

2.338 

1.01023 

1.529 

200 

1.01650 

3.248 

1.01096 

2.178 

was  found  to  be  satisfactory  for  the  range  of  f”  between  0  and  10.  For 
values  of  f"  greatei  than  10,  it  was  found  to  be  necessary  to  use  a  much 
smaller  spacing,  =  0.001,  cloae  to  the  wall. 


In  calculating  the  cases  for  both  compressible  and  incompressible  flows, 
it  is  necessary  to  transform  the  x,  y-coordinates  to  t,  rj-coordinates.  For 
flows  with  pressure  gradients,  it  is  also  necessary  to  calculate  the  pressure- 
gradient  term  P,  defined  by  (6.53).  The  5-  and  rj-coordinates  are  given  by 
(6.56)  and  (6.57),  respectively,  and  reduce  to 


x 

-/■ 


p  u  u  dx 
e  e  e 


(7.1) 


T|  = 


P  U 

e  e 


(p/pe)dy 


(7.2) 


for  two-dimensional  flows.  The  du  /dt-term  in  P  is  calculated  from  the 

e 

given  velocity  distribution  ug  versus  |  by  using  a  3-point  Lagrange  deriva¬ 
tive  formula 
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(4u./d‘>n  ■  ViV  .+  V.+  A”+iu<w 

n-i  n  n+i 

(7.3) 

where 

An-1 

*n  “  *n+l 

■  l  w  V(‘n-r  w 

(7.4a) 

A 

n 

8,n  "  {nn"  *n-l 
*n-l^*n~  *n+l^ 

(7.4b) 

An+1 

In” 

(*n+l"  *n-l^5n+l”  V 

(7.4c) 

For  compressible  flows,  the  fluid  properties  at  the  edge  of  the  boundary 
layer,  which  are  inputs  to  the  computer  program,  are  determined  by  the  test 
conditions.  On  the  other  hand,  for  incompressible  flows,  the  fluid  properties 
are  determined  at  standard  temperature  and  pressure. 

7.1  Incompressible  Flows 

The  test  cases  for  which  comparisons  are  calculated  by  the  present  method 
consist  of  a  variety  of  flows  of  widely  different  character:  flat-plate, 
flat-plate  with  mass  transfer  (blowing  and  suction),  equilibrium  flows  in  both 
favorable  and  adverse  pressure  gradients,  nonequilibrium  flows  in  adverse 
pressure  gradients,  separating  boundary  layers,  and  recovering  boundary  layers. 

7.1.1  Flat-plate  flows. 

Flat-plate  flow  offers  an  excellent  opportunity  to  compare  the  calculated 

results  with  experimental  results  as  well  as  with  well-established  correlations 

such  as  the  law  of  the  wall  and  the  defect  law.  A  flat-plate  flow  was  calcu- 

7 

lated  for  a  Reynolds  number  of  10  per  foot,  by  assuming  the  following  values 
of  the  flow  parameters: 

P  =  0.58  x  10-6lbf -sec/ft2 

P  =  0.237  x  10“2lbf-sec2/ft1+ 
ug=  1600  ft/sec. 
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The  value  of  u  was  selected  merely  as  a  matter  of  convenience.  For  this 

G 

case,  the  relations  given  by  (7.1)  and  (7.2)  reduce  to 


u  P  u  x 
e 

(7.5) 

P 

(7.6) 

respectively.  From  a  previous  study,  [7],  the  boundary-layer  thickness  was 

estimated,  and  the  approximate  value  of  at  each  (-station  was  determined. 

By  calculating  {  and  n  from  (7.5)  and  (7.6),  respectively,  a  flat-plate 

9 

flow  was  calculated  up  to  a  Reynolds  number  of  1  x  10  •  At  station  1,  at 

k 

Re  =  5  x  10  ,  the  flow  was  specified  to  be  turbulent.  The  calculated  value 
x  5 

of  Rg  at  Re^  =  1.05  x  10  was  identical  to  the  corresponding  value  of  R^ 
in  Table  II  of  [28].  The  initial  spacing,  h..,  was  0.001  and  K  was  1.01487* 

■L  g 

Even  though  the  calculations  showed  that  for  Rex>  5  *  10  the  valve  of 
should  iiave  been  about  150,  the  maximum  value  of  was  taken  to  be  100, 
because  of  the  computer-capacity  limit.  Approximately  70  stations  in  the 
(-direction  were  used.  The  time  per  station  was  approximately  10  seconds. 

The  results  of  the  calculations  are  presented  in  Figures  12a  through  12j. 

Figure  12a  shows  a  comparison  of  local  skin-friction  coefficients  calcu¬ 
lated  by  the  present  method  and  those  calculated  by  the  Prandtl-Schlichting 
formula  [29], 

cf  =(2  loe10Bex  -  0.65)"2*5  (7.7) 

as  well  as  experimental  values  and  Coles*  line  [28].  The  experimented  values 
are  taken  from  [28],  The  agreement  between  Coles’  line  and  experiment  is  very 

Q 

good.  At  very  high  Reynolds  number,  about  7  x  10  ,  the  calculated  values  of 
begin  to  deviate  slightly  from  Coles*  line.  This  is  believed  to  be  due  to 
the  smaller  value  of  ratner  than  150,  a  smaller  value  of  t^,  namely  100, 
was  used,  in  order  not  to  exceed  the  computer  capacity.  Figure  12b  shows  a 
comparison  of  calculated  shape  -factor  values  with  experimental  and  with  those 
given  by  Coles*  line.  In  figures  12c  and  12d  the  velocity  profiles  are  plotted 
in  law-of -the -wall  coordinates.  Figure  12c  shows  a  comparison  of  calculated 


results  with  experimental  and  with  those  given  by  the  universal  logarithmic 
velocity  # 

-h  -  5.75  log10  +  5.10 

U  AXf  V 

at  Re  ■  10^  and  10  .  Figure  I2d  shows  a  comparison  of  calculated  results 
with  those  given  by  Coles'  line.  In  figures  12e  and  12f  the  velocity  profiles 
are  plotted  in  velocity-defect  coordinates.  Figure  12e  shows  a  comparison  of 
calculated  results  with  experimental  and  with  those  given  by  the  logarithmic 
velocity  distribution[28] 


ue~u 

# 

u 


•80  -  5-75  log10C..05  -S-  ) 

o  e 


(7.9) 
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at  Re^  =  10  and  10'.  Figure  12f  shows  a  comparison  of  calculated  results 
with  those  given  by  Coles'  line.  Figure  12g  shows  a  comparison  of  calculated 
and  experimental  sublayer  profiles.  The  experimental  data  are  lCLebanoff's 
data  at  Re^  =  77,000  taken  from  [28].  The  agreement  of  results  in  figures 
12c,  12e,  and  12g  with  experiment  is  remarkable!  in  that  the  calculated  values 
agree  with  the  experimental  values  point  by  point. 

The  argument  of  the  logarithmic  term  in  (7*9)  is  sometimes  written  in 

4ta 

terms  of  y/A,  since  from  (6.123)  and  from  the  definition  of  u 


(7.10) 


It  can  also  be  written  in  terms  of  y/6.  Use  of  the  latter  presents  the 

difficulty  of  dealing  with  an  ill -defined  outer  edge  of  the  boundary  layer. 

On  the  other  hand,  ^6  =  3.6  [30]  for  constant-pressure  profiles.  For  flat- 

plate  flows,  6  is  therefore  related  to  the  defect-displacement  thickness  A 

by  o  A 

6  =  3^ 

* 

and  can  be  written  in  tersm  of  c„  and  6  as 


5  =  0.392837 


(7.11) 


v 
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Figures  12h  and  121  show  comparisons  of  calculated  and  experimental  [19] 

4 

shear-stress  distributions  and  mean-velocity  distributions  at  Re.  ■  7*7  x  10  , 
u  /u  ■  0.037*  The  definition  of  6  given  by  (7«ll)  is  used  for  the  boundary- 
layer  thickness  in  plotting  the  calculated  results.  The  shear-stress  coeffi¬ 
cients  given  by  (6.111*)  and  (6.115)  reduce  to 

c  «=  (2/e)1/2  u  f" 
tL 

c  =  (2/e)1/2  c  t" 

tT 

respectively.  The  agreement  of  the  velocity  profiles  is  nearly  perfect.  The 
shear-stress  distribution  is  also  very  good  except  near  the  edge.  This  is 
probably  due  to  the  ill-defined  definition  of  8  in  (6.26),  which  needs  to  be 
reconsidered.  For  example,  8  should  perhaps  be  based  on  the  point  where 
u*»  0.995*  Figure  12j  shows  a  comparison  of  calculated  eddy-viscosity  distri¬ 
butions  across  a  boundary  layer  with  values  calculated  from  both  Klebanoff's 
and  Townsend* s  data.  The  calculated  experimental  values  are  taken  from  [31]. 
The  results  indicate  that  the  calculated  eddy-viscosity  values  are  in  close 
agreement  in  the  inner  region  and  are  in  poor  agreement  in  the  outer  region. 


(7.12a) 

(7.12b) 


These  results  for  a  flat  plate  demonstrate  that  the  computing  program 
has  the  capacity  to  handle  Reynolds  numbers  greater  than  those  occurring  on 
any  existing  or  contemplated  aircraft,  as  well  as  on  most  ships. 


7.1.2  Flat-plate  flow  with  mass  transfer. 

The  present  method  is  readily  applied  to  flows  with  surface  blowing  or 
suction.  Mass  transfer  at  the  wall  is  handled  by  the  boundary  condition 
(6.59a),  which  reduces  to  / 


\(t  -  S0> 


n  ue(2£) 
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i  <  L 


i  >  t 


(7.15) 


for  uniform,  incompressible,  two-dimensional  flow.  In  (7.15),  £Q  represents 
the  distance  from  the  leading  edge  without  mass  transfer. 


( 


-64- 


7. 1.2.1  Uhlform  blowing  (Mickley  and  Davie).  The  experimental  data  used  for 
comparison  with  a  calculated  flow  with  uniform  blowing  were  obtained  by  Mlckley 

and  Davis  [32].  Two  blowing  rates,  v /u  -  0.001  and  0.002,  were  studied. 

™  e  = 

A  flat-plate  flow  was  calculated  for  a  Reynolds  number  of  3.123  x  10  per  ft 


by  using  the  values  of  p  and  p  that  were  used  in  Section  7.1.1>  with 
u  ■  50  ft/sec.  Thus,  an  effective  length  that  matched  the  momentum  thickness 
at  the  station  where  blowing  began  was  determined.  On  the  basis  of  this  infor¬ 
mation,  (now  lQ  is  known),  f  was  calculated. 


The  results  are  shown  in  figures  13a,  13b,  and  13c.  Figure  13a  shows  a 
comparison  of  calculated  and  experimental  momentum  thickness  and  local  skin- 
friction  coefficient  for  each  blowing  rate.  The  other  two  figures  show  com¬ 
parisons  of  calculated  and  experimental  velocity  profiles.  Figure  13b  shows 
the  calculated  and  experimental  values  for  y^/u  =  0.001  at  x  =  38.42  inches 
and  83*55  inches,  and  figure  13c  shows  tlie  calculated  and  experimental  values 

for  v  /u  «  0.002  at  the  same  stations.  The  agreement  between  the  calculated 
W'  e 

results  and  experimental  data  as  shown  by  the  figures  is  very  good. 


7»1. 2. 2  Uhifoim  suction  (Tennekes).  Experimental  data  obtained  by  Tennekes, 
[33]  and  [34],  were  used  for  comparison  with  a  calculated  flow  with  uniform 
suction.  Two  suction  rates,  v  /u  =  -  0.00312  and  —  0.00429,  were  studied. 

we  c 

A  flat  plate  was  calculated  for  a  Reynolds  number  of  8.13  x  Kr  per  ft  by 

using  the  values  of  p  and  p  that  were  used  before,  with  ug  =  130  ft/sec. 

Again,  as  with  blowing,  an  effective  length  was  determined  that  matched  the 

momentum  thickness  at  the  point  where  suction  began.  On  the  basis  of  this 

information,  f .  was  calculated, 
w 


Figure  l4a  shows  the  calculated  values  of  local  skin-friction  coefficient 


and  momentum  thickness,  together  with  the  experimental  data  deduced  from  the 
momentum  equation,  for  both  suction  rates.  Figure  l4b  shows  the  velocity-pro- 


file  comparisons  for  vw/ug  =  -  0.00312  at  x  =  1.25  feet  and  2.55  feet,  and 
figure  l4c  shows  the  velocity-profile  comparisons  for  =  -  0.00429  at 

x  =  0.594  feet  and  2.575  feet.  The  agreement  between  calculations  and  experi¬ 


ment  is  very  good. 
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7.1.3  Equilibrium  flows  in  both  favorable  find  adverse  pressure  gradients. 

All  the  previous  examples  were  constant  pressure  flows.  The  applicability 
of  any  general  method  depends  on  the  accuracy  of  the  results  it  gives  for  a 
wide  variety  of  flow  conditions;  hence  it  must  be  tested  on  flows  with  pressure 
gradients.  Indeed,  this  is  a  prerequisite  for  the  important  problem  of  pre¬ 
dicting  turbulent -boundary-layer  separation.  Because  of  the  interest  in  equi¬ 
librium  flows  that  are  characterized  by 

E  =  — ^  =  constant 

calculations  were  made  for  several  equilibrium  f J  ows .  Note  that  for  two-di¬ 
mensional,  incompressible  flows  with  pressure  gradients  the  (-coordinate  is 
determined  from  (7»l)>  which  reduces  to 

x 

(  =  UP J  ue  dx* 

0 

and  the  ^-coordinate  is  the  same  as  (7*  6).  The  following  sections  give  com¬ 
parisons  of  calculated  and  experimental  results  for  such  flows. 

7.1. 3.1  Equilibrium  flow  in  a  favorable  pressure  gradient  (Herring  end 
Norbury ) .  Herring  and  Norbury  [35]  experimentally  investigated  two  equilibrium 
flows  characterized  by  E's  of  —  0.35  and  -  0.53*  The  former  flow  was  calcu¬ 
lated  by  the  present  method.  A  free-stream  velocity  at  the  entrrnce  tc  the 

5 

test  section  of  72  ft/sec  and  a  Reynolds  number  of  4,5  x  l(r  were  thr  flow 
conditions  used.  Figure  15  shows  the  velocity  distribution  (obtained  by  trial 
and  error)  used  to  match  the  momentum  thickness  at  the  Initial  point  of  the 
experimental  data,  namely,  at  x  =  2  feet.  Note  that  with  this  effective 
length  added,  the  experimental  lengths  are  translated  4.37  feet. 

Figure  l6a  shows  a  comparison  of  calculated  and  experimental  momentum 
thickness  and  shape-factor  parameters,  together  with  the  experimental  velocity 
distribution,  and  figures  l6b  and  l6c  show  comparisons  of  cal'ruJateci  and  ex¬ 
perimental  velocity  profiles  at  x  =  2  and  3  feet  and  4  and  5  !cct;  respec¬ 
tively  .  The  agreement  with  experiment  is  very  good,  and  the  calculate*!  values 
of  velocity  profiles  correspond  to  experimental  values  poiut  by  point. 
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Figure  15.-  Velocity  distribution  usod  to  match  tho  momontum  thickness  for  the  experimental  data  of 
Herring  and  Norbury. 

7. 1.3*2  Equilibrium  flow  in  an  adverse  pressure  gradient  (Clauser).  In 
reference  [30],  Clauser  experimentally  obtained  two  equilibrium  flows  corre¬ 
sponding  to  two  different  pressure  distributions,  sometimes  designated  as 
P.D.l  and  P.D.2,  each  having  a  nearly  constant  value  of  E.  The  former  flow 
was  calculated  by  the  present  method.  A  free -stream  velocity  of  28.3  ft/sec 
at  the  beginning  of  the  flow  and  a  Reynolds  number  of  1.8  x  10  per  ft  were 
the  flow  conditions  used.  Figure  17  shows  the  velocity  distribution  used  to 
match  the  defect-shape  factor,  G,  of  the  experimental  data,  which  is  reported 
to  be  10.3*  The  defect-shape  factor  was  used  rather  than  the  momentum  thick¬ 
ness,  9,  because  the  momentum  thickness  was. not  reported  explicitly  in  [30]. 

Figure  18a  shows  a  comparison  of  calculated  and  experimental  shape -factor 
and  local  skin-friction  coefficient  parameters,  together  with  the  experimental 
velocity  distribution  taken  from  faired  velocity  distributions  given  in  [350]. 
Figure  l8b  shows  a  comparison  of  calculated  and  experimental  velocity  profiles, 
and  figure  l8c  shows  a  comparison  of  calculated  and  experimental  velocity 
profiles  in  the  defect-law  coordinates,  both  at  x  =  200  inches  and  375  inches. 
The  agreement  is  remarkably  good. 
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Figure  17.- Velocity  distributions  used  to  match  the  defect-shape  factor  for  the  experimental  data  of 
Clauser’s  P.D.1. 


and  Ferriss).  An  important  class  of  boundary  layer  flows  comprises  those  that 

have  been  perturbed  in  some  manner  and  are  then  allowed  to  recover  to  some 

equilibrium  condition.  As  an  example  of  this  kind  of  flow,  the  boundary 

layer  measured  by  Bradshaw  and  Ferriss[36]  was  considered.  A  portion  of  this 

-0.255 

flow  has  an  external  velocity  distribution  given  by  u^otx  *  from  x  =  0 
to  x  =  h'J  inches  and  is  followed  by  a  zero  pressure  gradient  from  60  inches 
to  x  =  95  inches.  A  reference  velocity  of  110  ft/sec  and  a  Reynolds  number  of 
6.85  x  10  per  ft  were  the  flow  conditions  used.  Figure  19  shows  the  velocity 
distribution  used  to  match  the  boundary-layer  parameters  at  x  =  23  inches. 

Note  that  when  an  effective  length  is  added,  the  experimental  lengths  are 
translated  by  2  feet. 

Figure  20a  shows  comparisons  of  calculated  and  experimental  momentum- 
thickness,  shape -factor,  and  local  skin-friction  parameters,  together  with  the 
experimental  velocity  distribution.  Figure  20b  shows  a  comparison  of  calcu¬ 
lated  and  experimental  velocity  profiles  at  47,  71  and  95  inches.  Even  though 
the  agreement  in  momentum  thickness  is  excellent,  except  for  the  first  station 
(x  =  23  inches )\  the  results  are  not  as  good  as  in  the  other  flows.  This  is 
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Figure  19.- Velocity  distribution  used  to  match  the  momentum  thickness  for  the  experimental  data  of 
Bradshaw  and  Ferriss. 

probably  due  to  the  inaccuracy  of  matching  the  initial  boundary-layer  para¬ 
meters  at  x  =  25  inches.  Generally,  several  trial-and-error  runs  (sometimes 
considerably  more)  are  necessary  to  match  the  initial  boundary-layer  parameters, 
which  are  usually  momentum  thickness  and  shape  factor.  At  least  one  accurate 
matching  of  one  of  these  parameters,  for  example,  momentum  thickness,  is  neces¬ 
sary  before  one  can  check  the  calculated  results  with  experiment.  This  is 
especially  true  in  an  equilibrium  boundary  layer  for  which  the  shape  factor  is 
nearly  constant  and,  unlike  the  shape  factor  of  a  nonequilibrium  flow,  it  does 
not  adjust  itself.  With  the  velocity  distribution  shown  in  figure  19,  the 
calculated  and  experimental  values  are 
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0.200 
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There  is  also  the  possibility  that  the  disagreement  was  caused  by  the  calcula¬ 
tion  of  P,  which  was  computed  by  the  formulas  given  in  Section  7.0.  Since  P 
involves  the  derivatives  of  velocity,  which  are  inputs  to  the  computer  program, 
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any  slight  irregularities  in  velocity  distribution  cause  inaccuracies  in  the 
0-values.  In  general,  whereas  some  small  irregularities  in  the  calculations 
of  0  do  not  appreciably  alter  the  results,  for  flows  that  are  near  separation, 
as  this  one  is,  irregularities  in  0  are  highly  undesirable.  For  this  reason, 
it  is  necessary  to  calculate  the  values  of  0  and  then  fair  the  results.  In 
order  to  improve  the  results  further  calculations  should  be  made. 

7.1. 4  Nonequilibrium  and  separating  flows. 

From  a  practical  standpoint,  nonequilibrium  and  separating  flows  are 
perhaps  the  most  important  flows,  since  they  are  often  encountered  in  the 
design  of  diffusers  and  lifting  surfaces.  For  this  reason,  three  separate 
flows  with  turbulent  separation  were  considered.  The  following  sections  give 
comparisons  of  these  flows  with  experimenta]  flows. 

7.1.^.!  Favorable  and  adverse  pressure  gradients  on  an  airfoil-like  body 
(Schubauer  and  Klebanoff).  In  reference  [37j>  Scnubauer  and  Klebanoff  experi¬ 
mentally  obtained  a  flow  characterized  by  an  initial  favorable  gradient  followed 
by  an  adverse  gradient  and  separation.  The  body  is  two-dimensional  and  has  a 
sharp  nose.  It  is  at  a  slight  angle  of  attack,  which  produces  a  pressure  peak 
at  the  leading  edge  that  causes  transition.  Separation  is  reported  to  have 
taken  place  at  25*7  feet  from  the  leading  edge.  A  reference  velocity  of 
loO  ft/sec  (velocity  at  17*5  ft  from  the  leading  edge)  and  a  Reynolds  number 
of  1  x  10  ^  per  ft  were  the  flow  conditions  used.  A  flat-plate  flow  was  used 
to  match  the  momentum  thickness  of  the  experimental  data  at  x  =  1  foot. 

The  results  are  shown  in  f  if  jure  s  21a,  21b,  and  21c.  Figure  21a  shows  the 
experimental  velocity  distribution  and  comparison  of  calculated  and  experimen¬ 
tal  momentum-thickness,  shape -factor  and  local  skin-friction  parameters.  The 
other  two  figures  show  comparisons  of  calculated  and  experimental  velocity 
profiles.  Figure  21b  shows  the  calculated  and  experimental  values  at  x  -  ^ 
feet  (in  favorable  pressure  gradient)  and  at  x  =  14.5  feet  (nearly  constant- 
pressure  flow).  Figure  21c  shows  the  calculated  and  experimental  values  at 
x  =  20  and  22  feet  (both  in  adverse-pressure-gradient  flow  and  conditions 
close  to  separation).  The  calculated  parameters,  with  the  exception  of  skin- 
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friction  ,  are  in  excellent  agreement  with  the  experimental  data  except  near 
separation.  The  present  method  predicts  separation  at  x  =  25  feet.  (Flow  is 
said  to  be  separated  at  the  last  station  before  the  one  at  which  f”  becomes 
negative;  when  f^  becomes  negative,  the  calculations  are  stopped.) 

7.1. ^.2  Adverse  pressure  gradient  on  a  body  of  revolution  (Moses).  Another 
example  of  a  nonequilibrium  boundary  layer  with  adverse  pressure  gradient  is 
a  flow  studied  by  Moses  [38],  which  was  used  as  a  basis  of  comparison.  The 
pressure  distribution  designated  P.D.2  was  considered.  The  results  reported 
are  for  the  cylindrical  section  only.  Because  the  body  studied  is  a  cylinder 
of  constant  diameter,  it  is  not  a  test  of  the  method  for  axisymmetric  flow, 
since  with  the  first  theory  and  with  the  present  transformations,  the  flow  is 
two-dimensional.  The  experimental  separation  point  was  at  about  28.8  inches. 

A  free -stream  velocity  of  35*3  ft/sec  at  the  beginning  of  the  flow  and  a 

5 

Reynolds  number  of  2.2  x  10  per  ft  were  the  flow  conditions  used.  Figure  22 
shows  the  velocity  distribution  (obtained  by  trial  and  error)  used  to  match 
the  momentum  thickness  at  the  initial  point  of  the  experimental  data,  namely, 
x  =  0. 


Figure  22.- Velocity  distribution  used  to  match  the  momentum  thickness  for  the  experimental  data  of 
Moses'  P.D.2. 

* 

Personal  communication  with  Klebanoff  indicated  the  experimental  values  of 
these  quantities  to  be  in  error. 
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Note  that  when  an  effective  length  is  added,  the  experimental  lengths  are 
translated  by  9»1  feet. 


Figure  23a  shows  comparisons  of  calculated  and  experimental  momentum- 
thickness,  shape-factor,  and  local  skin-friction  parameters,  together  with 
the  experimental  velocity  distribution.  Also  shown  in  this  figure  are  the 
two  separation  points  determined  by  the  present  method.  The  calculated  re¬ 
sults  are  based  on  two  values  of  7  in  the  convergence  criterion  for  f^, 
namely, 


.pll  _  x>»l 

1Q+1  ZQ 


<  7 


The  results  indicated  by  the  solid  line  in  figure  23a  are  those  that  are  ob¬ 
tained  by  specifying  7  to  be  0.01,  and  the  results  indicated  by  the  dotted 
line  are  those  that  are  obtained  by  specifying  7  to  be  0.C01.  Except  near 
the  separation  point,  the  calculated  results  are  identical  and  are  in  excellent 
agreement  with  the  experimental  values.  According  to  the  present  method,  sepa¬ 
ration  based  on  7  ~  0.001  is  indicated  at  x  =  27*5  inches,  and  separation 
based  on  7  =  0.01  is  indicated  at  x  =  29  inches.  Figure  23b  shows  a  compari¬ 
son  of  calculated  (7  =  0.001)  and  experimental  velocity  profiles  at  11,  20, 

26,  and  29  inches.  Again,  except  near  separation,  the  calculated  results  are 
in  excellent  agreement  with  experiment.  Figure  23c  shows  a  comparison  of 
calculated  (7  =  0.001)  and  experimental  turbulent  shear-stress  distributions 
at  11,  20,  and  2h  inches.  The  agreement  is  fair. 


7.1.  *+.3  Adverse  pressure  gradient  on  an  airfoil  (Von  Doenhoff  and  Tetervin), 

A  further  example  of  a  nonequilibrium  boundary  layer  with  adverse  pressure 
gradient  is  the  NACA  65(216)-  222  airfoil  tested  by  von  Doenhoff  and  Tetervin 
[39]*  The  pressure  distribution  at  an  angle  of  attack  of  10.1°  and  a  chord 
Reynolds  number  of  2.6^  x  10^  was  considered.  The  experimental  separation 
point  was  at  about  55-percent  chord.  A  free-stream  velocity  of  211  ft/sec  at  the 
beginning  of  the  flow  and  a  Reynolds  number  of  1.315  x  10^  per  ft  were  the  flow 
conditions  used.  A  flat-plate  flow  was  used  to  match  the  momentum  thickness 
of  the  experimental  data  at  7. 5"P^rcent  chord,  rather  than  the  leading  edge  as 
the  origin.  This  was  necessary,  since  the  boundary  layer  was  tripped  at  7‘5- 
percent  chord. 
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The  results  are  shown  in  figures  24a,  24b,  24c,  and  24d.  Figure  24a 
shows  the  experimental  velocity  distribution  and  comparisons  of  calculated 
and  experimental  momentum-thickness,  shape-factor,  and  local  skin-friction 
parameters.  The  other  three  figures  show  comparisons  of  calculated  and  ex¬ 
perimental  velocity  profiles  up  to  the  separation  point,  namely,  at  55-percent 
chord.  In  figure  24a,  the  calculated  momentum  thickness  values  are  in  good 
agreement  with  the  experimental  values.  However,  the  agreement  is  poor  for 
the  shape-factor  values.  This  is  somewhat  surprising,  because  the  calculated 
velocity  profiles  are  in  good  agreement  with  experiment  (for  example,  at  25- 
percent  chord).  Since  the  shape  factor  is  the  ratio  of  displacement  thickness 
to  momentum  thickness,  a  better  correlation  in  H-values  would  normally  be 
expected.  The  calculated  velocity  profiles  begin  to  deviate  considerably  from 
the  experimental  values  near  separation.  Figure  24d  shows  the  calculated 
velocity  profiles  at  50-  and  54-percent  chord  with  the  separation  profile  at 
55-percent  chord.  The  present  method  predicts  separation  at  54. 5-perc.ent 
chord,  which  is  in  excellent  agreement  with  the  experimental  separation  point. 

7.2  Compressible  Flows 

The  test  cases  for  which  comparisons  are  calculated  by  the  present  method 
consist  of  flat-plate  flows  with  and  without  heat  transfer.  Only  a  few  test 
cases  were  calculated,  mostly  because  of  the  difficulties  experienced  in 
obtaining  convergence  of  the  numerical  solution  of  the  momentum  equation.  As 
a  result,  most  of  the  time  was  devoted  to  attempts  to  stop  the  oscillations  and 
obtain  convergence.  In  fact,  at  present  the  numerical  method  for  compressible 
flow  still  needs  to  be  refined.  The  calculations  for  compressible  flow  use 
5-point  averages  of  the  coefficients  defined  by  (6.78)  in  the  rj-direction. 

This  scheme  was  initially  used  to  stop  the  oscillations,  and  because  of  the 
limited  time  the  computer  program  was  not  refined.  That  is,  the  averaging 
process  in  the  T)-direction  was  not  confined  only  to  the  eddy-viscosity  terms, 
as  it  was  in  the  incompressible  flows. 

7.2.1  Flat-plate  flow  with  adiabatic  wall. 

Three  separate  sets  of  experimental  data  were  used  to  test  the  present 
method  for  a  fxat-plate  flow  with  adiabatic  wall  conditions.  The  first  set  of 
experimental  data  are  those  of  Spivack  [4o]  at  Me  =  2.8.  The  measurements 
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were  made  on  the  flat  wall  of  a  two-dimensional  nozzle.  The  flow  was  calcu- 

7 

lated  for  a  Reynolds  number  of  2  x  10  per  ft,  by  assuming  the  following 
values  of  the  flow  parameters: 

U  =  0.376  x  10~^lb„-sec/ft2 
e  I  1 

p  =  0.23^  x  10”2lb_-sec2/ft^ 
e  f  ' 

u  =  3170  ft/sec 
e  ' 

Te  =  550  deg  R 
H  =  8.2  x  106ft2/sec2 

The  value  of  u  was  calculated  from  the  definition  of  M, 
e  ’ 

u  =  49.1  M 
e  e  *  e 

and  the  total  enthalpy,  Hg,  was  calculated  from 


For  this  case,  the  transformed  x-  and  y-coordinates  are  given  by  the  relations 
(7.1)  and  (7.2),  respectively. 


Figures  25a  and  25b  show  the  calculated  local  skin-friction  parameters  and 
experimented,  values  that  were  obtained  from  the  momentum  equation,  namely, 


dx  2 


The  parameters  are  plotted  against  Reynolds  number  based  on  length,  Rex,  in 
figure  25a  and  on  Reynolds  number  based  on  momentum  thickness,  Re^,  in  figure 
25b.  The  agreement  is  very  good  in  both  cases. 


The  second  set  of  experimental  data  are  those  of  Chapman  and  Kester  [hi] 
at  Me  =  0.8l,  2.5,  and  3.6.  The  experimental  setup  consisted  of  a  cylindrical 
model  with  air  flowing  axially  at  the  outer  surface  of  the  model.  The  experi¬ 
mental  skin-friction  values  were  obtained  by  a  direct  measurement  of  toteil 
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force  on  the  cylinder.  The  flow  was  calculated  by  assuming  the  following 
values  of  the  flow  parameters: 

6  P 

Ue  =  0.376  x  10“  lb^-sec/ft 

pg  =  0.234  x  10“2lbf-sec2/ft^ 

Tg  =  530  deg  R 

The  values  of  u  and  H  for  each  Mach  number  were  calculated  in  the  same 
e  e 

way  as  in  the  preceding  case. 

Figure  26  shows  a  comparison  of  calculated  and  experimental  average  skin- 
friction  coefficients  at  the  three  Mach  numbers.  The  agreement  is  excellent 
for  Mg=  0.8l,  but  only  fair  at  Me  =  2.5  and  3.6. 

The  third  set  of  experimental  data  are  those  of  Matting  et  al  [42]  at 
Mg  =  0.20,  2.95 1  and  4.2.  The  measurements  were  made  on  the  flat  wall  of  a 
two-dimensional  nozzle.  The  calculations  used  the  values  of  the  flow  parameters 
tabulated  below. 


M 

e 

6 

p  x  10 
e 

lbf-sec 

ft2 

p  x  102 

He 

..  2 
lbf-sec 

~~nr~ 

T 

e 

deg  R 

0.20 

0. 376 

0.234 

530 

2-95 

0.17° 

0.55 

219 

h.2 

0.110 

0.93 

132 

The  values  of  u  and  H  for  each  Mach  number  were  calculated  in  the  same 
e  e 

way  as  in  the  preceding  cases. 

Figure  27  shows  comparisons  of  calculated  and  experimental  local  skin- 
friction  coefficients  at  the  three  Mach  numbers.  The  agreement  for  Mg  =  0.20 
is  excellent,  but  lor  the  other  two  Mach  numbers  is  only  fair.  Comparisons  of 
calculated  and  experimental  boundary-layer  Mach  number  profiles  and  velocity 
profiles  for  Me  =  2.95  and  Rex  =  31  x  10^  are  shown  in  figure  28a  and  for 
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M  =4.2  and  Re  =  69  x  10  in  Figure  28b.  The  agreement  in  velocity  profiles 

0  X 

is  better  than  the  agreement  in  Mach  number  profiles.  However,  in  both  cases 
the  agreement  is  only  fair. 

7.2.2  Flat-plate  flow  with  heat  transfer. 

Only  one  set  of  experimental  data,  those  of  Pappas  [43],  was  used  for  a 
flat-plate  flow  with  heat  transfer.  The  test  Mach  number  used  in  the  compari¬ 
son  is  I.69.  The  experimental  skin-friction  values  were  obtained  from  the 
momentum  equation.  Two  different  ratios  of  wall  temperature  to  edge  tempera¬ 
ture,  namely  Tw/Tg  or  gw,  were  considered.  The  calculations  used  the  follow¬ 
ing  values  of  the  flow  parameters: 

p  =  0.376  x  10  ^lb„-sec/ft2 
e  f  ' 

p  =  0.234  x  lC_2lb„-sec2/fti+ 
e  f  ' 

Tg  =  550  deg  R 

The  values  of  ug  and  Hg  were  calculated  in  the  same  way  as  in  Section  7«2.1. 

Figure  29  shows  comparisons  of  calculated  and  experimental  average  skin- 
friction  coefficients  for  gw  =  1.70  and  1.6l.  Figure  30  shows  comparisons 
of  calculated  and  experimental  local  Stanton  num^rs  for  the  same  g^values. 

The  agreement  in  both  cases  is  fair.  Figure  31  shows  a  comparison  of  skin- 
friction  variation  with  Mach  number.  The  ratios  of  local  compressible  skin- 
friction  values  to  incompressible  values,  namely  cf/cf  >  were  calculated  at 
several  Mach  numbers  and  compared  to  some  experimental  results.  The  agreement 
seems  to  be  satisfactory. 
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8.0  CONCLUDING  STATEMENTS 

A  numerical  solution  of  the  turbulent -boundary  layer  equations  based  on 
an  eddy-viscosity  concept  and  an  eddy-conductivity  concept  are  presented.  The 
numerical  method  used  to  solve  the  boundary  layer  equations  differs  from  the 
one  used  in  a  previous  study  [7]  under  contract  NOw  64-0352c.  In  reference 
[7],  the  validity  of  using  an  eddy-viscosity  concept  in  solving  the  equations 
of  the  incompressible  turbulent  boundary  layer  was  explored.  The  momentum 
equation  was  solved  in  its  nonlinear  form  by  the  old  shooting  method  of  pre¬ 
vious  studies  of  laminar  flows  [8  through  13].  The;  shooting  method  had  proved 
itself  in  accuracy  and  reliability  for  laminar  flows,  and  it  was  therefore 
preferred  to  other  numerical  methods,  which  are  usually  based  on  linearized 
solutions.  The  previous  study  [7]  indicated  that  the  eddy -viscosity  concept, 
when  used  with  an  accurate  technique  for  solving  the  boundary-layer  equations, 
was  capable  of  predicting  results  that  were  in  good  agreement  with  experiment 
and  many  of  the  characteristic  features  of  the  turbulent  boundary  Jayer.  The 
main  o'  ’ections  to  that  method  were  long  computing  times  and  restrictions  on 
the  step  sizes  in  the  streamwise  direction.  In  the  former  case,  for  example, 
a  typical  test  case  consisting  of  about  23  stations  in  the  x-direction  took 
approximately  one  hour  on  the  IBM  7094  computer.  That  amounted  to  about  2 
minutes  per  station.  The  numerical  method  used  in  that  study  also  had  the 
disadvantage  that  very  short  steps  in  the  streamwise  direction  —  which  are 
essential  near  separation  —  could  not  be  taken. 

For  these  reasons,  when  the  efforts  were  extended  to  the  solution  of  the 
equations  of  the  compressible  turbulent  boundary  layer,  it  was  necessary  to 
use  a  different  numerical  method.  Almost  the  same  eddy-viscosity  formulation 
as  in  [7]  was  used,  except  for  small  modifications  tc  account  for  the  compressi¬ 
bility  effects.  The  results  indicate  that  the  present  method  is  very  fast.  A 
typical  computation  time  for  incompressible  flows  on  the  IBM  709^  computer  is 
approximately  10  seconds  per  station,  which  means  that  the  present  method  is  at 
least  10  times  as  fast  as  the  old  method.  A  typical  computation  time  for  com¬ 
pressible  flows,  on  the  other  hand,  is  about  15  seconds  per  station.  In  addi¬ 
tion  to  having  very  short  computing  times,  the  present  method  has  no  restric¬ 
tions  on  the  step  sizes  in  the  streamwise  direction;  that  is,  step  sizes  can 
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be  as  small  as  necessary. 

Various  incompressible  and  compressible  turbulent  flows  are  calculated 
by  the  present  method,  and  comparisons  with  experiments  are  made.  Use  results 
for  incompressible  flows  are,  in  general,  very  good,  and  so  far  the  results 
for  compressible  flow  are  encouraging.  More  test  cases  need  to  be  run  for 
both  flows,  especially  for  compressible  turbulent  flows.  The  results  presen¬ 
ted  do  not  represent  a  finished  development,  but  are  only  what  has  been  ob¬ 
tained  by  one  particular  formulation  of  eddy  viscosity  and  turbulent  Prandtl 
number. 
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9.0  TECHNOLOGICAL  FORECAST* 

9.1  Problem 

A  method  of  calculating  velocity  and  temperature  profiles,  skin  friction, 
displacement  thickness,  momentum  thickness,  heat  transfer,  and  the  separation 
point  of  the  turbulent  boundary  layers  is  presented.  Accurate  calculation  of 
these  quantities  can  be  extremely  beneficial  to  the  design  of  ships,  turbo¬ 
machinery,  heat  exchangers,  lifting  surfaces,  and  all  types  of  aerospace 
vehicles. 


9.2  State  of  the  Art,  Solution,  and  Forecast 

[See  Section  5»0  of  this  report]. 

9.3  Suggestions  and  Implications 

The  present  method  uses  eddy-viscosity  and  eddy-conductivity  concepts  in 
obtaining  a  numerical  solution  of  the  turbulent -boundary-layer  equations.  The 
method  of  solution  has  the  advantage  of  being  applicable  to  both  laminar  and 
turbulent  flows.  Since  most  flows  have  a  laminar  portion  in  the  vicinity  of 
a  stagnation  point  or  leading  edge  (for  example,  on  an  airfoil  at  high  Reynolds 
numbers),  such  an  advantage  is  particularly  useful.  The  present  method  uti¬ 
lizes  the  upstream  history  as  the  calculations  proceed  downstream.  At  any 
station,  a  laminar  profile  can  be  obtained  by  merely  setting  the  eddy-viscosity 
term  in  (6. 5l)  equal  to  zero. 

The  present  method  also  has  the  advantage  that  various  formulations  of 
eddy  viscosity  or  turbulent  Prandtl  number  can  be  used  with  very  little  change 
in  the  basic  method.  Since  the  results  presented  in  this  report  are  obtained 
by  one  particular  formulation  of  eddy  viscosity  and  turbulent  Prandtl  number, 
such  an  advantage  may  be  very  useful  in  improving  the  results  for  both  incom¬ 
pressible  —  where  the  results  are  very  good  except  near  separation  —  and  com- 
* 

This  section  of  the  report  is  included  in  response  to  the  requirements 
of  the  contract. 
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pressible  flows.  For  example,  in  the  outer  region,  the  present  method  uses 
an  eddy-viscosity  relation  modified  by  an  intermittency  factor  y  to  describe 
the  turbulent  shear  transport.  This  relation,  suggested  by  Clauser  [5]  for 
equilibrium  boundary  layers,  is  given  by  (6.27) 

eo  =  p  k2Ue5*7 

where  k2  is  a  constant  of  proportionality.  For  equilibrium  boundary  layers 

this  constant  is  0.018  [5].  In  both  incompressible  and  compressible  studies 

here, however,  this  constant  was  taken  to  be  0.0168,  the  value  given  in  [18], 

In  this  way,  (6.27)  was  also  extended  to  nonequilibrium  boundary  layers.  In 

order  to  improve  the  results  near  separation,  it  is  quite  possible  that  the 

* 

velocity  and  length  scales,  namely  ug  and  6  ,  in  (6.27)  may  need  to  be 
changed.  According  to  a  recent  study  [44],  it  appears  that  when  different 
length  and  velocity  scales  are  used  in  normalizing  the  eddy-viscosity  formula 
for  the  outer  region,  somewhat  less  variation  in  eddy  viscosity  from  station 
to  station  is  obtained  than  when  (6.27)  is  used.  The  results  of  these  calcu¬ 
lations  for  the  experimental  data  in  [44]  are  summarized  in  Table  III,  which 
presents  the  maximum  and  minimum  values  for  four  different  normalizations  of 
eddy  viscosity.  The  first  is  used  in  this  report  and  evidently  is  the  worst 
of  the  lot. 

TABLE  III.  EDDY  VISCOSITY  VARIATIONS 


e 

e 

e 

e 

* 

u  6 

*  * 
u  6 

u  6 

* 
u  e 

e 

e 

Maximum  value 

0.028 

0.79 

0.059 

1.4 

Minimum  value 

0.0048 

0.26 

0.014 

0.57 

Max/Min 

5.85 

3.04 

2.78 

2.75 

Another  modification  of  the  outer  eddy-viscosity  relation,  which  may  im¬ 
prove  the  results,  is  the  use  of  a  mixing  length  concept,  as  in  the  inner 
region.  Escudier  and  Spalding  [4^]  have  recently  studied  the  distribution  of 
mixing  length  in  the  turbulent  boundary  layers.  The  results  indicate  that  for 
the  outer  region  of  the  boundary  layer  the  mixing  length  is  constant  and  is 
equal  to  0.075  &• 
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The  results  for  compressible  flows  indicate  that  it  is  also  necessary 
to  make  improvements  in  the  formulation  of  the  turbulent  Prandtl  number.  Even 
though  the  velocity  profiles  agree  quite  well  with  experiment,  the  temperature 
profiles  show  considerable  deviation  from  experimental  values  for  compressible 
flows,  especially  close  to  the  wall.  These  results  suggest  that  it  is  neces¬ 
sary  to  use  a  separate  expression  for  turbulent  Prandtl  number  in  each  region 
and  that  the  initial  assumption  of  a  constant  turbulent  Prandtl  number  is  not 
very  satisfactory,  since  the  variation  of  turbulent  Prandtl  number,  especially 
close  to  the  wall,  is  considerable.  Figure  32  shows  turbulent  Prandtl  numbers 
of  boundary  layers  on  a  cooled  plate  [46].  This  figure  supports  the  necessity 
of  expressing  turbulent  Prandtl  number  separately  in  each  region. 
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Figure  32.* Turbulent  Prandtl  numbers  of  boundary  layers  on  a  cooled  flat  plate,  Ma^~  5.1,  (Te-Tw)/  Te'x  0.35. 


In  summary,  it  can  be 
couraging,  the  work  should 
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said  that  because  the  results  have  been  so  en- 
be  continued  along  these  lines. 
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appendix  A 


BEHAVIOR  OF  THE.  EDDY  VISCOSITY  VERY  CLOSE  TO  THE  WALL 

The  inner  region  of  a  turbulent  boundary  layer  contains  a  layer,  commonly 
called  the  sublayer,  adjacent  to  the  wall,  where  the  flow  is  primarily  viscous, 
since  all  velocities,  including  turbulence  fluctuations,  become  zero  at  the  wall. 
This  region  is  not  uniform  either  with  respect  to  time  or  with  respect  to  dis¬ 
tance  along  the  wall.  In  the  region  y  >  5^,  where  6^  is  the  sublayer  thick¬ 
ness,  the  effect  of  viscosity  in  the  flow  decreases  gradually  with  increasing  dis¬ 
tance  from  the  wall  until  it  finally  becomes  negligibly  small.  Beyond  this 
point  the  flow  is  completely  turbulent.  The  intermediate  region,  where  the  flow 
is  neither  completely  viscous  nor  completely  turbulent,  is  called  the  transition 
region  [31]. 


In  describing  the  behavior  of  the  eddy  viscosity  very  close  to  the  wall, 

4 

one  can  use  the  above  model  and  can  show  that  e  should  vary  as  y  as  y  -»  0 
provided  disturbances  in  the  boun  lary  layer  are  two-dimensional  and  periodic. 
Since 


—  u'v* 


£  du 
P  5y 


(Al) 


and  if  bn/by  is  assumed  to  be  constant  in  the  viscous  sublayer  (this  corre¬ 
sponds  to  a  linear  velocity  distribution),  it  is  only  necessary  to  show  that 
—  u*v’,  which  is  proportional  to  e,  should  vary  as  y  as  y  -♦  0. 


Consider  a  two-dimensional  incompressible  mean  flow  and  a  likewise  two- 
dimensional  disturbance.  Assume  that  the  mean  velocity  U  depends  only  on  y, 
that  is,  U  =  U(y),  as  in  parallel  flows.  It  is  also  necessary  to  assume  that 
the  pressure  in  the  main  flow  depends  on  x  as  well  as  on  y,  because  the 
pressure  gradient  bP/bx  maintains  the  flow.  It  can  be  shown  that  under  these 

* 

The  flow  in  the  boundary  layer  can  also  be  regarded  as  a  good  approximation 
of  parallel  flow,  because  the  dependence  of  the  velocity  U  in  the  main  flow 
on  the  x-coordinate  is  very  much  smaller  than  the  dependence  on  y  [29  ]• 
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conditions  the  Navier-Stokes  equations  and  the  continuity  equation  can  be 
written  as  [29] 


du* 

5 T 


+  u 


du1 

s r 


+  v* 


1  V 
p  5x 


(A2) 


dv' 

ST 


+  u 


i 

p 


(A3) 


du*  ,  dv* 

sr  f 


0 


(A4) 


where  the  primes  denote  the  fluctuations.  By  differentiatirc  (A2)  with  re¬ 
spect  to  y  and  (A3)  with  respect  to  x  and  subtracting  the  two  equations, 
the  fluctuating  pressure  term,  p',  can  be  eliminated.  The  resulting  expression 
can  be  written  as 


St-  usr 


(A5) 


where  o>  is  the  vorticity  component  perpendicular  to  the  plane  of  motion, 
z 

It  is  defined  by 


dv*  du* 

“z  =  5T  “  5T 

Next  introduce  a  stream  function  representing  a  disturbance  of  the  form 

+(x,  y,  t)  =  <p(y)  eia^X  "  Ct^  (A7.) 

where  <p  is  the  amplitude,  a  is  a  real  quantity  defined  by  L,  A  is 
the  wavelength,  L  is  a  characteristic  length,  and  c  is  a  complex  number  de¬ 
fined  by  a  propa3ation  velocity  c^  and  an  amplification  factor  c^,  that  is, 

c  =  c  +  i  c . 

r  l 


Now,  very  close  to  the 
may  be  writter  as: 


wall  (y  -»  0),  where  U 


v 


and 


v*  are  very  small,  eq.  (A5) 


(A8) 
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Since 


»•  =  |  =  <p'(y)eia(x-ct> 

(A9) 

v*  *  -  |  =  -  i  «  9(y)ela(x-ct) 

(A10) 

the  vorticity  component  can  be  written  as 

m  =  c?  _  T(y)e«»(^rt)  (An) 

z 


by  using  (A 6).  Substituting  the  expression  for  u)  in  (All)  into  (A8)  and 

z 

simplifying  gives 


2 
d  T 

dy2 

The  solution  of  (A12)  is 


i  Cite 


T  =  0 


(A12) 


T  =  c^e 


i  '\liotcr/v'  y  +  c  e-i  Viacr/v*  y 


which  reduces  to 


T  =  c^e 


i  Viator/v* 


(A13) 


if  T  is  assumed  to  be  finite  when  y  -»oo.  Equation  (A15)  can  be  expanded  in 
the  form 


T  =  c. 


V"3c  1  Oc 

2^r  y  “  i  2T  y  +  •  •  * 


(Al4) 


Substituting  the  expression  (Al4)  into  (All)  and  neglecting  the  terms  higher 


than  y  in  (Al4)  gives 

y"  -  Q29‘?'  cjl  -  (l  -  i)V-2V  y  “  i 


Ocr  21 

"2V  yJ 


(A15) 


The  solution  of  (A15)  is,  using  the  method  of  variation  of  parameters, 


qp  =  A  e 


^  +  Be-^y  _  4 


_  /Oc  '  Oc  0  _ 

1  -  (1  -  ijy  -5“  y  -  i  “2^y2  + 


(A16) 
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The  constants  A  and  B  can  be  evaluated  from  the  boundary  conditions,  namely, 
v*  =  u*  =  0  at  the  wall.  Therefore  it  follows  from  (A9)  and  (A10)  that  <?*  (y) 
and  <p(y)  are  equal  to  zero  at  the  wall.  If  the  constants  A  and  B  are 
evaluated  in  this  way,  (Al6)  now  becomes 

<p  =  y  y2  -  5  U  -  A  (a2  -  2i  Jy14  (A17) 

The  fluctuating  velocities  u*  and  v*  can  now  be  obtained  by  using  (A9)  and 
(A10)  and  by  considering  the  real  part  of  the  resulting  expressions,  that  is. 


u'  =  Re  (d*/dy)  ,  v*  •■=  — ■  Re  (d*/dx) 


(A19) 


Therefore,  by  using  (Al8)  and  (A19)  it  can  be  shown  that 


and  that  should  vary  as  y11  as  y  -*  0,  where  n  £ 
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APPENDIX  B 

MODIFICATION  OF  THE  INNER  EDDY-VISCOSITY  EQUATION  FOR  COMPRESSIBLE  F. 

The  expression  for  the  inner  eddy  viscosity  given  by  (6.23)  can  also  be 
used  for  compressible  flows  if  p  is  taken  to  be  a  variable  and  if  the  ex¬ 
ponential  term  is  modified  to  account  for  the  heat  transfer  in  the  sublayer. 

A  logical  generalization  is  to  consider  a  Stokes-type  flow  with  the  following 
conditions : 

1.  Sinusoidally  oscillating  infinite  wall. 

2.  Variable  p,  T,  p  . 

3.  Oscillations  small  enough  so  that  any  compressibility  effects  they 
cause  can  be  neglected. 


Because  of  the  last  condition,  the  energy  equation  does  not  enter  into 
the  problem.  For  a  nonsteady  parallel  flow,  the  Navier-Stokes  equation  re¬ 
duces  to  (see  eq. (5.17)  of  [29]) 


du 

St 


a  / 
57 


du 

¥ 


) 


Introducing  the  transformation 


(Bl) 


dy  =  p  dz 

(B2) 

into  (Bl)  yields 

chi  d2u 

(B3) 

To  solve  (B3),  let 

U  =  fM  ol(ot  *  c) 

(bM 

Then 

i  of  el(w-*  c) 

(B5) 

d2u  i(ot  ♦  c)  d2f 

Sz2  =  5  fc2 

(B6) 

Substituting  the  expressions  given  by  (B5)  and  (B6)  into  (B3)  and  rearranging 
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gives 


-p  |i  i  of  =  0 


Assume  that  pp  is  averaged  over  a  certain  distance,  and  write  the  average  as 
pp  .  Then  the  solution  of  (B7)  can  be  written 

f  =  A  exp  (l  +  i  )V“h'  z  +  B  exp  -a  + (B8) 

since  Vi™  =  <• 

V2 

Assume  that  f  remains  finite  as  z  -*  oo  so  that  A  =  0.  Then  substitu¬ 
ting  the  resulting  expression  in  (B^)  into  (Bb)  yields 


u  =  B  exp  — 


PP  0 


+  i(ot  +  e  — 


■VvVi 


(B9a) 


Then  (B9a)  can  be  written  as 


u  =  B  exp(— 


z)  cos  ,  at  +  c  — 


(B9b) 


or  by  using  the  transformation  given  by  (BP),  eq.(B9t)  becomes 


u  =  B  exp  (-V  V  J  * 


COS  ot 


* '  -Vv/  f } 
0 


(bye) 


The  expression  given  by  (B9>')  can  be  approximated  by  expanding  p  by  a  Taylor 
series  around  an  average  value  of  e,  that  is, 


=  P  +  y(de/d y)+  . . . 


(BIO) 


If  the  viscosity  in  (BIO;)  is  introduced  into  (Bye)  and  if  the  se  -ond-order 
and  higher  order  terms  of  y  are  neglected,  (Boo)  can  be  written  as 


u  -  B  ex 


COS  O  t  +  € 


-W- 


(sii) 
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Equation  (fill)  is  the  same  as  the  incompressible  solution  except  that  mean 
values  of  p  and  p  are  used.  A  mean  value  looks  more  accurate  then  the 
wall  value.  The  mean  value  of  v  should  be  that  corresponding  to  the  sublayer 
thickness. 


Since  the  cosine  term  is  bounded, 


(Bll)  can  be  written  as 


(B12) 


The  purpose  of  modifying  the  mixing  length  and  consequently  the  expression 
for  eddy  viscosity  given  by  Prandtl's  formula,  (6.21),  is  to  account  for  the 
.  ublayer  and  to  account  for  the  transition  region  between  the  sublayer  and  the 
fully  turbulent  part  of  the  inner  region.  For  this  reason,  Van  Driest  [lU] 
defines  the  mixing  length  as 


/  =  k^ytl  -  exp(—  y/A)] 


Because  the  Stokes  effect  is  localized  to  a  region  very  close  to  the  wall, 
v  and  p  should  assume  values  appropriate  to  this  region.  A  convenient  ex¬ 
pression  can  bo  formed  by  introducing  wall  values,  but  corrected  by  a  factor 
(v„/7)^2  to  account  for  the  small  deviation  from  wall  values,  as  follows: 


exp, 


-  Vt? 


=  exp  1  - 


_o 

r2v 


V 


exp 


-VF 


\ 


(B13) 


The  mixing  length  now  becomes 


1  —  exp 


i-VFi) 


t  =  ^  y 
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APPENDIX  c 

VARIABLE -GRID  DIFFERENTIATION  FORMUIAS 


Hie  first-,  second-,  and  third-derivative  formulas  for  the  variable -grid 
system  discussed  in  Section  6. 6. 1,3  can  be  obtained  in  the  following  way: 
Consider  the  Lagrange  interpolation  formula 

n 

f(x)  =  ^  ii(x)f(xi)  +  E(x)  (Cl) 

i=0 


where 

^j(x) 


t(x) 

(x  -  xi) 

(x  -  xq)  ...(x  -  x._1)(x  -  xi+1)  ...  (x  -  xn) 
(x.  —  X  ) . . .  (x X  .  ,)(x.—  X.,,)  . . .  (x.  —  x  ) 

i  o'  l  i-l,v  l  l+l '  '  i  n' 


and 


f(»+l) 


(0 


Differentiating  (Cl)  gives 

n 

fr(x)  /[(x)  f(x.)  +  Er(x)  (C2) 

i=0 

The  equations  (Cl)  and  (C?)  are  now  evaluated  for  three-  and  five-point  systems 
to  obtain  the  first -derivative  formulas  for  the  variable -grid  system  with  three 
points  and  the  first-,  second-  and  third  derivative  formulas  for  the  variable- 
grid  system  with  five  points.  Figures  Cl  and  C2  show  the  variable-grid  systems 
for  three  and  for  five  points,  respe  -lively. 


First-Derivative  Formulas  for  Three  Points 


fj-l(x)  "  4.[-  B X(P  +  Kjf.^-  E?(I  +  K)f.  -  B  f.+1]  1/.  f3(|)  (©) 


K  1m  , 
I 


n  M  =  h.r-  B1K  f  +  Bg(l  -  h)  f.  +  B,fi+1)  f'(0 


(CM 


Figure  C2  *  Variabla-grid  system  for  five  points. 
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f'+1(x)  =  h.tB^  timl+  B2(l  +  K)f.+  B5(l  +  2K)f.+1]  +  2(0  (C5) 


where 


B,  =  - 

1  h*  (l  +  K) 


(C6a) 


b2  =  --T" 
hf  K 

l 


(C6b) 


B,  =  - 

^  K(l  +  K) 


(C6c) 


First  Derivative  Formulas  for  Five  Points 
Definitions  of  new  terms  follow  the  formulas 


fj.2(x) 


=  —  h?  +  a.aT+  a_az+  a.a^az)  f.  *  +  B_a_a0azf.  . 

1-2  112  15  25  12  5  i-2  2  1  2  5  l-l 


B3S2a3fi  +  Bljaia’.Ii*l  *  Bsaia2ri*'^  +  12c  ^  *'i ^ 


hi-2^  B1K  ala21i-2  +  a2+  ala2  K  ala?^fi-l 

a  VVj':  +  VViu*  B/aifi+?]  I'*) 


fj  (x) 


hi.2K\tBiK?‘'i.2+  V  aifi-it  k  V  V 


.  * 


B'*aifi+1  B,/i+P-*  +  IPO  l‘i-2  1  ^ 


fi+i(x)  =  hi-rK^“  BiKVi-2  -  Wi-i  -  V  v/i  +  MaiVaiarK 


-a/-aiK  )fi+1  +  - 


K°a  aQ  i, 

h._g  HO  (cio) 


I 


f?+2(x)  -  hi_2K  [B;La1a2K  fi_2+  B2aiajK2fi_1+  33a2a3K  fi+  \aia2a3fi+l 

K^8>  &  GL 

+  B^Ca^K  +  a^K^  a^K  +  ^^3)^+2 3  +  120  ^  hi-2f5^ 


(cn) 


Second-Derivative  Formulas  for  Five  Points 


fi-2(x)  =  2hi-2[Bl(ai+  V  V  aiV  aia3+  a2a3)fi-2+  B2(aia2+  aiV  a2a3)fi-l 

+  B5(a2+  a^+  a2a5)f.+  B,t(a1+  a^+  a1a3)fi+1+  B^a^  a2+  aia^'fi-i-2^ 

^ala2+  aia3+  a2a3+  ala2a3^hi-2f (C12) 


fj.^x)  =  ^_2K[B1K(a1+  a2+  B2(aia2K  +  agK  +  a^K  +  ap+  a^f^ 

+  B3(Kaia2"  a2"  ai)fi+  W  a2’  l)fi+l+  B5(Kar  V  l)fi+2] 

+  K2 (a^+  a2+  a^-  t)  (C13) 


q  (x) 


=  2h2_2K[B1K2(Ka1-  a^-  l)f._2+  B^OC2-  a^-  1)^+  B^-  ^ 

-  K2-  K2a1+  K^f.f  B^d  -  Kax-  K2 )i*i+1+  B^-  Ka^  K2 )fi+2] 

-  £  K'a^l  -  a  -  Kax+  K^^U)  (ClU) 


fj+1(x)  =  2h2_2K[B1K2(a1-  Ka^  K2)^*  B^K(a?-  Ka,-  K3)^*  B^a^- 
-  K5ai)f.+  B^a^  Ka-  K^-  K^-  K*)fi+1+  B^a^*- 

+  K2 a1)f.+2  +  ^  K3(aia2-  a^K  -  a^-  ajK3)^^ |)  (C15) 
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fV+2(x)  =  2hJ_2K[B1K2(a1a2+  Ka2+  *?\)fi+2+  B^^ay-  Ka^  K^f.^ 

+  B5(a2a3+  ^  a3+  Va2a3+  Kala5+  ^^i+i 

+  B5(a2a?+  Kaxa?+  l^a^f  K2 a^a^  K3a2+ 

+  K^(aia2K^+  aia3K2+  a2a3K  +  aia2a3^hi  (Cl6) 

Third-Derivative  Formulas  for  Five  Points 


fjl'(x)  = 


-  6hi-2[Bl(l  +  V  V  a3)fi-2+  B2(ai+  a2+  a3)fi-l 


+  B5(l  +  a2+  a3)f.+  B^l  +  a1+  a5)fi+1+  B^l  +  &1+  a2)f.+2] 

+  ^(ax+  a2+  a?+  axa2+  a2a5)hj  2  f5(|)  (C17) 


f*'»(x)  =  -  ^.^[B^d  +  &1+  a2)f.  _2+  B2(Ka2+  Ka^ 

+  B3(Ka2+  Kax-  l)f.+  B^Ka^  K  -  l)f.+1+  B^Ka^  K  -  l)f.+2] 

+  ^  KCKa^  Ka2+  Ka^  a?+  a^h^f^O  (Cl8) 

f}*'(x)  =  6  fei_2tB1K(l  -  K  -  Ka^f.^*  B2(ai-  K2- 

+  B^a^  K  -  K2-  A1)fi+  Bu(&1+  K  -  1^)1^+  B5(ax+  K  -  K2)^] 
+  ^  K(a]—  Ka^-  Ka2-  K- K^h2^  f5(0  (C19) 


-99- 


f}”(x)  =  6h.  ^[B^a^  K  -  K^)fi_2+  B2(a2+  K2-  K5)^  B3(a2+  K^-  K3)^ 

+  B4(a2+  Ka1+  K2-  K3)f.+1+  B^a^  Ka^  K2)^] 

+  ]“  K(ajLa2+  Ka2-  K^-  K3a^—  KU)h2_2  f5(ft)  (C20) 

f”*(x)  =  6h._2[B1K(a2+  Ka1+  K2)^^  B2(a3+  K2&1+  K3)^^  B^a^  Ka2+  K3)f± 
+  B4(a5+  Ka2+  K^Jf^  +  B^a^  Ka2+  K3)f.+2] 

+  ^  K(a2a3+  Kaxa3+  K2  a,  +  K^a^  K3a2+  Ka^  h2_2  f5(ft)  (C2l) 

The  constants  B^»  B2,  By  3^,  B,.,  and  a^,  a2,  &y  are  f^iven  by 


Bi  =  T - 

(C22a) 

hi-?ala?a3 

(C22b) 

(C22c) 

(C22d) 

B*  *  T“T - 

"l-s'W 

(C2?e) 

1  +  K  +  K2,  a3  =  l  +  K+  K2+K3 

(C22f ) 

ax  =  1  +  K, 


Figure  12.- Results,  of  calculations  for  a  flat-plate  flow. 

(a)  Comparison  of  calculated  local  skin-friction  coefficient  with  experiment.  Prandtl-Schlichting  formula,  and  Coles'  line 
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?  igur®  12. -Continued,  (c)  Comparison  of  calculated  and  experimental  law-of-the-wall  profiles. 


Figure  12.- Continued,  (d)  Comparison  of  calculated  law-of-the-wall  profiles  with  Coles*  line. 


Figure  12.- Continued,  (e)  Comparison  of  calculated  and  experimental  velocity-defect-law  profiles. 


Figure  12.- Continued,  (f)  Comparison  of  calculated  velocity-defect-law  profile*  with  Cole*'  lino. 


Figure  12.- Continued,  (g)  Comparison  of  calculated  and  experimental  sub-layer  profile  according  to  Klebanoff. 


Figure  12. -Continued,  (i)  Comparison  of  calculated  and  experimental  mean-velocity  distributions  according  to  Klebanoff. 

Re,  =  7.7«104.  u* /u,  =  0.037. 
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CALCULATION,  REF.  31 
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Rag  s  7.7  x  104,  ~~  =  0.037 
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Figure  12.- Continued,  (j)  Comparison  of  calculated  eddy-viscosity  distributions  across  a  boundary  layer 

with  values  calculated  from  Klebanoffs  and  Townsend's  data. 


Flgur»  13-  Continued,  (b)  Comparison  of  calculated  and  experimental  valocity  profile*  for  v*/u#  =  0.001  at  stations  x  =  38.42  and  83.56  inches. 


Figure  13.- Continued,  (c)  Comparison  of  calculated  and  experimental  velocity  profiles  for  vw/ue  =  0.002  at  stations  x  =  38.42  and  83.55  inches 
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Figure  14.- Results  uf  calculations  for  tha  flat  plate  with  two  uniform  suction  rates,  vw/ut  =  -0.00312  and 
-0.00429  (experimental  data  of  Tennekes). 

(a)  Comparison  of  calculated  and  experimental  momentum  thickness  and  local  sk  n-f fiction 
coefficient  for  the  two  suction  rates. 
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Figure  14. -Continued,  (b)  Comparison  of  calculated  and  exoerimental  velocity  profiles  for  vw/ue  = -0.00312  at  stations  x  =  1.25  and  2.55  feet. 


x  ~  FEET 


Figure  16. -Results  ot  calculations  for  an  equilibrium  flow  in  a  favorable  pressure  gradient  with  E  =  -0.35 
(Experimental  data  of  Herring  and  Norbury). 

fj)  Experimental  velocity  distribution  and  comparison  of  calculated  and  experimental  momentum 
thickness  and  shape-factor  parameters. 
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Continued,  (b)  Comparison  of  calculated  and  exper. menial  velocity  profiles  at  stations  x  =  2  and  3  feet. 


Figure  16.-Continued.  (c)  Comparison  of  calculated  and  experimental  velocity  profiles  at  stations  x  =  4  feet  and  5  feet 
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Figure  18.- Results  of  calculations  of  an  equilibrium  flow  in  an  adverse  pressure  gradient  (Experimental 
data  of  Clauser's  P.D.1). 

(a)  Experimental  velocity  distribution  and  comparisons  of  calculated  and  experimental 
shape-factor  and  local  skin-friction  parameters. 


Figure  18. -Continued,  (c)  Comparison  of  calculated  and  experimental  velocity  profiles  in  the  defect-law 

coordinates  at  x  =  200  and  375  inches. 


Figu*  20.- Continued,  (b)  Comparison  of  ca  leu  la  tad  and  axparimantal  valocity  prof  I  las  at  x  =  47,  71  and  95  inchss. 
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Figure  21. -Results  of  calculations  for  favorable  and  adverse  pressure  gradients  on  an  airfoil-like  body 
(Experimental  data  of  Schubauer  and  Klebanoff). 

(a)  Experimental  velocity  distribution  and  comparisons  of  calculated  and  experimental 
momentum-thickness,  shape-factor,  and  local  skin-friction  parameters. 
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Igure  25.- Comparison  of  calculated  and  experimental  local  skin-friction  coefficients  for  an  adiabatic 
(Experimental  data  of  Spivack). 

(a)  Variation  of  local  skin-friction  with  Reynolds  number  based  on  the  length,  Rox 


(Experimental  data  of  Chapman  and  Kester). 


PRESENT 


Comparison  of  calculated  and  experimental  local  skin-friction  coefficients  for  an  adiabatic  flat  plate  at  M  =  0.20,  2.95  and  4.20 
(Experimental  data  of  Matting  et  al). 


Figure  29.- Comparison  of  calculated  and  experimental  average  skin-friction  coefficients  for  a  flat  plate  with  heat  transfer  at  M  =  1.69 
(Experimental  data  of  Pappas). 
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